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The form factor K{t) is calculated analytically to the order as well as numerically for a 
rectangular billiard perturbed by a 5- like scatterer with an angle independent diffraction constant, 
D. The cases where the scatterer is at the center and at a typical position in the billiard are 
^— ^ . studied. The analytical calculations are performed in the semiclassical approximation combined 

' with the geometrical theory of diffraction. Non diagonal contributions are crucial and are therefore 

^SJ , taken into account. The numerical calculations are performed for a self adjoint extension of a 5 

I > ■ function potential. We calculate the angle dependent diffraction constant for an arbitrary perturbing 

O ' potential U{r), that is large in a finite but small region (compared to the wavelength of the particles 

, that in turn is small compared to the size of the billiard). The relation to the idealized model of 

the (5-like scatterer is formulated. The angle dependent diffraction constant is used for the analytic 
^ ^ ' calculation of the form factor to the order r^. If the scatterer is at a typical position, the form factor 

^ . is found to reduce (in this order) to the one found for angle independent diffraction. If the scatterer 

' is at the center, the large degeneracy in the lengths of the orbits involved leads to an additional 

r \ \ small contribution to the form factor, resulting of the angle dependence of the diffraction constant. 

■ The robustness of the results is discussed. 

Ch ■ I. INTRODUCTION 

(N ■ 

^ I The spectral properties of quantum mechanical systems are closely related to the ones of their 

\Q . classical counterparts [0,0- This relation can be made precise and quantitative in the semiclassical 
O ' 

' limit. In the extreme cases this correspondence is quite well understood, although some important 



O 
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> 



details are subject of research The spectral statistics exhibit a high degree of universality, 

although there are exceptional systems that do not exhibit this universal behavior, due to some 
special properties. When a classical system is chaotic the statistical properties of corresponding 
quantum spectrum are well described by Random Matrix Theory (RMT) with the symmetry 



^ . determined by the underlying Hamiltonian. For integrable systems, on the other hand, levels are 
d ' typically uncorrelated resulting in Poissonian spectral statistics . There are however intermediate 
situations. One is of mixed systems, where in some parts of phase space the motion is chaotic and 
in other parts it is regular leading to different statistics . Another type of systems that exhibit 
intermediate statistics are systems with singularities, that are integrable in absence of these singu- 
larities. These are models of physical situations where an integrable system is strongly perturbed 
in a region much smaller then the wavelength of the quantum particle. Examples of such systems 
are billiards with flux lines, sharp corners and 5-like interactions [0,^,^, n , T^,]!^ . The speciflc 



system that will be analyzed in the present work is a billiard with a 5-like potential perturbation. 
The interest in billiards of various types is primarily theoretical since it is relatively easy to analyze 
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them analytically and numerically. Billiards were studied also experimentally for electrons p^ , p^ 



microwaves [0,^,^ and for laser cooled atoms p2| , p3| . We hope that in the future, perturbations 
of the type discussed in the present work will be also introduced in the experimental realizations of 
billiards. 

Semiclassical methods provide a convenient path to explore how classical dynamical properties 
affect the quantum spectral statistics. The connection between spectral statistics and classical 
dynamics is usually made using trace formulas which connect the quantum density of states to a 
sum over classical periodic orbits. For classically chaotic systems the trace formula was developed by 
Gutzwiller |2^j2^ while for classically integrable systems the trace formula was developed by Berry 



and Tabor |2g]. Such formulas can be used to compute the energy level-level correlation function, 
R2{r])i where 77 is the level separation. Its Fourier transform is the spectral form factor K{t), where 
the "time" r is conjugate to rj. For systems with strong perturbations on length scales smaller 
than the wavelength, standard semiclassical theory cannot be used and diffraction effects have to 
be taken into account. For billiards with localized perturbations this can be done in the framework 
of the Geometrical Theory of Diffraction (GTD) p^. In this approximation, scattering from the 
perturbation can be described by a free propagation to the perturbation multiplied by a diffraction 
constant, which depends on the incoming and outgoing directions, and then free propagation away 
from the perturbation. 

The spectral statistics of the rectangular billiard with a flux line was recently studied |T^,|TB|. In 
this system the flux line affects the amplitudes of the contributions of periodic orbits. When the form 
factor is computed it is found that at small times r it tends to a flux dependent constant (instead of 
1 in absence of the flux). Certain triangular billiards were shown to have similar behavior The 
works mentioned above are based only on the contributions from periodic orbits. It is of interest 
to include the contributions of orbits that are diffracted by the flux line in the spectral form factor. 
The contribution of diffracting orbits turns out to be rather complicated. The contribution of orbits 



with one diffraction was computed by Sieber |2^, while Bogomolny, Pavloff and Schmit computed 
the contributions of diffracting orbits with several forward diffractions to the density of states [p9| . 
The effect of diffraction on the spectral statistics in these systems was not determined yet. The 
difficulty results from the need to know the contributions of orbits which diffract almost in the 
forward direction and to include all of these in the form factor, and this has not been done so far. 
Since the diffraction from flux lines and corners is complicated it is of special interest to examine 
a system in which the diffraction is relatively simple. A class of such systems, that consists of 
integrable billiards with localized perturbations, is the subject of the present work. For systems 
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of this type only diffraction effects are responsible for deviation of the behavior from the one of 
integrable systems, therefore these are most appropriate for the exploration of such effects. 

The diagonal approximation, where the cross terms between the contributions of classical orbits 
are ignored because of rapidly varying phases associated with them, is used for a large variety of 
systems. Berry has shown that using the diagonal approximation leads to agreement with the form 
factor of random matrices, at small times [^]. Only recently Sieber and Richter were able to 
go beyond the diagonal approximation and to compute the next power of r for chaotic systems with 
time reversal symmetry. This non diagonal contribution was found to originate from orbits with 
small angle self intersections and orbits which are close to them in coordinate space, except near 
the intersection. Such configurations are not possible if there is no time reversal symmetry, thus 
this is an example where the difference in the spectral statistics is connected to a difference in the 
dynamical properties of the classical counterpart, and the difference can be obtained directly from 
a formula involving the classical periodic orbits. For chaotic systems with singular perturbations, 
where effects of diffraction are important, also the non diagonal contributions to the form factor 
K{t) are important. The final result, however, is that such perturbations do not affect the spectral 



statistics as is clear from the works of Sieber p2| , |33| and of Bogomolny, Leboeuf and Schmit ||3j 

In contrast to the case of chaotic systems, addition of a localized scatterer, is expected to change 
qualitatively the spectral properties of integrable systems. Moreover, for such systems the diagonal 
approximation is not sufficient to obtain any interesting physics and one should go beyond this 
approximation. Non diagonal contributions can be calculated following a method that was developed 
by Bogomolny By this method contributions of some manifolds consisting of orbits that are 

almost parallel are calculated. If a point interaction, that is a generalization of a ^-function potential. 



is added to a rectangular billiard, one obtains the so called Seba billiard |[TT|] . The spectral statistics 
of this system were computed exactly, subject to the assumption that the underlying integrable 



system exhibits Poissonian statistics |]13[ and it was shown not to be those of integrable or chaotic 
systems. The spectral statistics of the Seba billiard with periodic boundary conditions were shown 



by Berkolaiko, Bogomolny and Keating to be identical to those of certain star graphs, where 
the form factor was calculated by Berkolaiko and Keating 



After the work was completed a paper by Bogomolny and Giraud |^, where the expansion of 
K{t), was calculated to all orders in r was published on the Web. It is more general than the present 
paper, but in our paper more of the physical properties are calculated and discussed in detail. The 
methodology of the derivation is different. A comparison between various results of these two papers 
is presented in App. E. 
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Generalized 5-functions were extensively studied in mathematical physics in the framework of 
the self adjoint extension theory . The calculations in |13|,|6| 



were 



performed within this framework. In the present work also a potential of finite small extension a is 
studied, under the assumption that a is much smaller than the wavelength of the quantum particle, 
that in turn is much smaller than the dimensions of the billiard. The relation of the (5-like scatterers 
to well defined potentials is discussed in the paper. 

The aim of this work is to use periodic and diffracting orbits to compute the form factor of rect- 
angular billiards with localized perturbations. The calculation includes non diagonal contributions 
and leads to a power expansion of K{t) where only the first three powers of r are calculated. The 
structure of this paper is as follows. In Sec. |I| the form factor is calculated for angle independent 
diffraction. In Sec. |IT| the model of point 5-like interactions is presented, in order to introduce a 
simple example of angle independent scattering. For this model the results of Sec. || are compared 
to numerical results in Sec. The scattering from localized potentials that are physically well 
defined is studied in Sec. 0, where it is shown how to compute the diffraction constant in the limit 



where the size of the perturbation is smaller than the wavelength. In Sec. ^ the form factor for 
angle dependent scattering, with the diffraction constant that was calculated Sec. 0, is computed. 
The results are summarized and discussed in Sec. IVU. 



II. THE FORM FACTOR OF A POINT-LIKE S WAVE SCATTERER 



In this section the form factor will be calculated for a rectangular billiard with a scatterer such 
that its diffraction constant D, is angle independent. The form factor is 

/oo 
drjR2{v)e''''^\ (1) 
"OO 

where R2{i]) is the energy levels correlation function: 

R2{v) = {dose [e - ^) d^sc [e + ^)) a' (2) 

and A denotes the mean level spacing of the system. The brackets denote averaging over an energy 
window which is large compared to the mean level spacing A but small when compared to the 
semiclassical energy E. In the semiclassical limit, the oscillatory part of the density of states, 
dosc{E) = d{E) - {d{E)), is of the form 

(iosc(^) ^ E +C-C., (3) 

p 

where c.c. denotes the complex conjugate and the sum is over periodic orbits. Using this density of 
states, the form factor is given by: 
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K{t) = 27r^A ApA;,e^(^''-V)5 I ^ 
\ pp' \ 



where Sp is the action of the orbit, tp is its period, and r > is assumed. The sum is over the 
leading semiclassical contributions to the density of states. For smooth potentials dose is given by 



the Berry- Tabor formula for integrable systems |26| and by the Gutzwiller trace formula for chaotic 



systems [p^pSf . In these formulas the sum is over the periodic orbits of the system. In the present 
work, periodic orbits and diffracting orbits will be included in the sum. 

Due to the averaging only action differences that are of the order of h or less can contribute to 
the form factor. As a result, the diagonal approximation can be used, where only pairs of orbits 
with the same action contribute. This approximation is valid for integrable systems and thus 
also valid when there are only once diffracting orbits (since their density is proportional to that of 
periodic orbits as will be explained in the following). When orbits with more than one diffraction 
are included this approximation is no longer valid and non diagonal terms should be included as is 
done in what follows. Since the motion in the billiard is free, the periods of the orbits are related 
to the lengths by tp = ^ where k = \fE is the wavenumber (the units = 1, 2m = 1 are used 
from now). The mean level spacing , in the semiclassical limit, is given by A = In the diagonal 
approximation the form factor is expressed as a sum over orbits of identical lengths 

Kir) = j:^{\Apmp-Akr)). (5) 

p 

The averaging widens the delta function and enables to replace the sum over orbits by an integral 
using the density of orbits. Every factor of length, Ip, is replaced by a factor of Akr, and as will be 
seen, every diffraction order eventually contributes a power of r to the from factor. The form factor 
K{t) will be calculated to the third order in r. 

The spectral statistics are very sensitive to the location of the scatterer, since they are influenced 
by degeneracies in the lengths of the diffracting orbits. Therefore we choose to treat first the simplest 
case when the scatterer is in the center of the rectangle. Then length degeneracies are maximal. 
Later we treat the case in which the scatterer is at a typical location, such that its coordinates 
(relative to the rectangle sides) are typical irrational numbers. 

A. Scatterer at the center 

The periodic and diffracting orbits are classified by the number of bounces from the boundaries. 
We will denote by Np = {Np, Mp) a periodic orbit with 2A'p and 2Mp bounces from the boundaries 
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and by Nj = {Nj,Mj) a diffracting orbit with Nj and Mj bounces from the boundary. When the 
scatterer is at the center of the rectangle the length of the diffracting segments is given by 



= ^alNf + alMf (6) 

while the length of the periodic orbits is 



l, = 2^alN^ + alM^ (7) 

where and (Xy cirG the lengths of the sides of the rectangle. 

The contributions to the oscillatory part of the density of states were calculated in the framework 
of geometrical theory of diffraction in App. A. For an angle independent scatterer the contribu- 
tions (|A25|) to dose can be simplified since D{6, 6') reduces to a constant D and if the scatterer is at 



the center the terms do not depend on ^ and u. The contribution of them with no more than three 
diffractions is 



where 



4°) = ^=e-f, (9) 



AY' = ^ ' ' De-'-^\ (10) 



and 



37rA;(87rA;)i Jl~l~i~, 



We omitted the factors fij defined in ([A14|) since the fraction of orbits where these differ from unity 



is negligible in the calculations of the present section. Examples of the various orbits contributing 
to (11-0) are presented in Fig. |l|. 
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FIG. 1. Several orbits contributing to the density of states: a) a periodic orbit (the (1, 1) orbit), b) a forward once diffracting 
orbit (2,2), c) a once diffracting orbit (2, 1), and d) a twice diffracting orbit. 



First the contribution of terms that involve no more than one diffraction, resulting from the first 
two terms in the expansion (P) of dose-, will be calculated. Following ( |A17|) the contributions of the 
diffracting terms can be classified according the parity of Nj and Mj. If both Nj and Mj are even 
(even-even orbits) the orbit j coincides with a periodic orbit, since in the notation of ( A17 ), at the 
scatterer 6 in = 9 out (compare Fig. |l]a) and b)) and the contribution of these orbits will be calculated 
first. There is no such relation for the other parities, (even-odd, odd-even, odd-odd, referring to 
the parities of Nj and Mj) and actually for such parities the diffracting orbits are halves of periodic 
orbits (this is only true if the scatterer is at the center, therefore for example the orbit in Fig. [I| c) 
is not a half of a (2, 1) periodic orbit). For even-even diffracting orbits the first two sums in the 
expression (H) can be grouped as 

d^ol{E) = Ape""'^ + ^f^''''' + c-c. (13) 
p j 

where the first sum runs over periodic orbits while the second sum runs over diffracting orbits, 
omitting orbits where both Nj and Mj are even and 



A, = I + ^!^ De-t I e-^t. (14) 
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The diagonal approximation is now apphed to various terms in (^) and the sum in (|^) is replaced 
by an integral with the density 

pW = ^ (15) 

of orbits of length /. This is the density of periodic orbits as well as the density of diffracting orbits 
with a given parity of Nj and Mj as can be easily seen from (^) and (^. The resulting contribution 
of the periodic and even-even diffracting orbits to the form factor is 

^(i,e)(^) = r dli^ + ^.^D + S{1 - Akr) 

= l + S3Dr + ^r2. (16) 

The contribution of the other diffracting orbits (odd-odd, odd-even, even-odd) with the amplitudes 
Af^ of (|10|) are 

i^(i.°-°)(r) = ir(i'°-^)(r) = K^^''-°\t) = (17) 
and their combined contribution is 

ir(^-)(r) = ^4^. (18) 
Thus, the contribution of the once-diffracting orbits and the periodic orbits to the form factor is 

ir«(r) = ir(i'^)(r) + ir(^'°)(r) = 1 + '^Dt + \D\W (19) 



From (p|-|Tl|) we note that the powers of and ^Jlj (the powers in (|9|-p!TD combined with the 
increasing density of multiple diffracting orbits resulting from powers of (|1^)) increase with the 
powers of D and by (^) and (^ these lead to increasing orders of r. We are interested first in all 
contributions to the order r^. The additional terms that may contribute to this order result from 
non diagonal terms in (^ where p is a periodic orbit and j' is an orbit with two diffractions, namely 
terms of the form 

ApA*,e^^^^-^^'^ = ApAf*/^^^^-^n-^^^\ (20) 

The contribution of such terms will be analyzed in the following. All other contributions are of 
higher order in r. 

We turn to analyze the contribution of non diagonal terms involving periodic and twice diffracting 
orbits. These oscillate as e^'^'^'^'^'^i^'^a) that may be potentially fast. Naively, one can guess that after 
smoothing, the contributions relevant to the form factor are those for which Ip = Ij^ + Ij^. However, 
it turns out that there are not enough orbits that satisfy this equality in order to compensate for 
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the extra factor of due to the additional diffraction. The contribution comes from combinations 
of orbits for which the lengths are almost equal. Such combinations survive since there are enough 
orbits for which the phase of the exponent is of the order of unity in the semiclassical limit. This is 
caused by the fact that these orbits have hidden saddle manifolds on which their action differences 
are slowly varying. These manifolds were identified by Bogomolny and were used to compute 
the density-density correlation function and mean values of powers of Green's functions for integrable 
systems. The diffraction problem discussed here has lengths (or actions) that are related to the length 
of the periodic orbits so that saddle manifolds exist. We use Bogomolny 's method to compute this 
contribution. It will be convenient to expand the length of the diffracting orbits around the length 
of orbits that satisfy Ip = Ij^ + ■ This expansion is given by 



+ , ^ ialN.5N. + alM.6M.) + rfM.^iV,- - (21) 

^alNf + alMf ' ' v ^ ^> 2{alN'j + alM^)i ^ ^ ^ ^> 

The length difference in the exponent is 

5lp = Ij, + Ij^ - Ip (22) 

The periodic orbit is not expanded because for each periodic orbit there will be a saddle manifold 
composed of different twice-diffracting orbits. The summation over the orbits in the manifold will 
be also combined with summation over all of the periodic orbits. Equivalently, it is possible to fix 
any of the lengths, expand the other lengths, and to reach the same contribution to the form factor. 

The condition that the orbits are on the saddle manifold is that the length difference, 6lp, is 
stationary leading to: 

1 . o , , . 1 



alNJN,, + a^M,,5M,,) + ^==^===(alNJN^, + alM^JM,,) = 0. (23) 



alNl + alMl " JalN^ + alMl 



The orbits on the manifold are the integer solutions of this equation. Since the aspect ratio of the 
rectangle is assumed to be a typical irrational number, such solutions can be found only if the orbits 
are related by 



alNl + alMl 

(24) 



alNl + alMl r,. 



This means that the orbits are repetitions of some simpler orbit, for example 



^alNl + alMl = r,,^aln^ + alm^ = r.Jo- (25) 
Substituting this condition in (|2^) leads to 
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alnSNj^ + alnSNj^ + almSMj^ + almSMj^ = 0. (26) 
Since the aspect ratio is a typical irrational number an integer solution can be found only if 



5% = (27) 

We are interested in the contribution that survives the energy averaging, required for contributions 
to the form factor. Therefore we will consider the saddle manifold for which 

I, = 2^alNl + a2M2 = (r,, + r,J/o. (28) 

This condition states that the periodic orbit (which will contribute) is also proportional to (n, m) 
and Ip = 2rpZo with 

2rp = rj, + rj, (29) 

is satisfied. 

The contribution of the terms of this saddle manifold to the level-level correlation function R2{v) 

is: 

+C.C. + (//-> -r]) . (30) 

The symbol [r] —rf) denotes two terms similar to the first ones but with the opposite sign of rj. 
This is actually a sum over {Np,Mp), {Nj^^Mj^) and {Nj.^,Mj^) that define the orbits p, ji and j2 
respectively. The contributions with similar phases, that do not cancel each other, and therefore 
survive the averaging over energy involved in the calculation of the form factor are those that 
satisfy (|27|) and (^). The sum will be calculated and then Fourier transformed in order to obtain 
the corresponding contribution to the form factor. The dominant contributions are these in which 
the length difference 6lp is much smaller than unity. The configuration of these contribution is shown 
in Fig. |. 
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FIG. 2. The configuration of lengtlis on tlie saddle manifold 



The two diffracting segments are almost parallel to the periodic orbit. Such a configuration is 
allowed since the difference between any lattice points is also a lattice point, moreover, the condi- 



tions (|27|) , (|29|) force the orbit to have this configuration. Note that /i,/2 are the projections of the 
lengths of the diffracting segments on the direction of the periodic orbit (to the leading order in 
S(f)j). Since they are almost parallel they are a good approximation to the lengths of the diffracting 
segments. 

Since the contributing orbits have the configuration of Fig. |^, the variables li and t are convenient 
coordinates. This choice also takes into account the conditions { ^7\ ) and (P^). The length difference 
is given by 

1 2 /I 1\ 1 4 /I 1' 



= + U + ilJ-g' U^lJ (31) 

where I2 = Ip — h- The length difference Sip appears as a phase in an exponent multiplied by k. Thus 
only length differences which scale clS CcLll contribute to the correlation function. The lengths of 
orbits that will contribute to the form factor scale as k and thus the corresponding values of t are of 
the order of unity or less. Since the sum is dominated by terms in which t ~ 1 or less, the expansion 
(pip is justified. This expansion may break down when either li or I2 ^ 0. However, there is a 
minimal length for the diffracting orbits and such lengths do not occur in the sum. 

The correlation function is computed by replacing the sum over orbits by an integral over li and 
t. The slowly varying terms are replaced by their values on the saddle (when t = 0) and the 
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phase is replaced by the leading term of (|31|). To transform the sum into an integral the density 
of contributions is needed. The configuration in Fig. ^ is obtained when the plane is tiled with 
rectangles using reflections with respect to the hard walls. The scatterer is also reflected and its 
images form a lattice (since there is a scatterer in the center of each rectangle). An example for this 
unfolding is presented in Fig. 13 (for a scatterer which is not at the center). 




FIG. 3. Unfolding of orbits using reflections around the rectangle boundaries. 



The (1, 2) orbit in the original rectangle (solid line) is equivalent to the one connecting the scatterer 
with its image (dashed line). The density of points in this lattice is -j. After these manipulations 
the correlation function is given by 



+ c.c + {r] ^ -7]). 



Calculation of the Gaussian integral leads to 



The Fourier transform leads to the form factor 



(32) 



(33) 



(34) 
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where only positive times were taken into account. The sum over periodic orbits is replaced by an 
integral using the density ([T5|) leading to 

K'^^\t) = T^^{-D^), (35) 

which is the non diagonal contribution to the form factor. This contribution is indeed of order as 
was expected. 

After computing all of the contributions of order we turn to compute the contributions to 
the form factor which are of order r^. There are three possible contributions of this order: diagonal 
contributions from twice diffracting orbits that will be denoted by K^^\t), non diagonal contributions 
resulting from combinations of once and twice diffracting orbits, denoted by K'^^'^^t), and from 
combinations of periodic orbits and orbits with three diffractions, denoted by K'^^''\t). These 
contributions are computed in the following. 

To obtain the diagonal contribution of twice diffracting orbits (||) is substituted in (H), and only 
the terms with two diffractions are kept. A Fourier transform leads to 

K('){r) = / A(2;A(2^e-^'^('^i+'--'«-'-)5 (bl±In±Iil±Il± _ ^\ \ 4nkA (36) 
\j1j2j3j4, \ ^ / / 

where contributions from negative times were omitted. The diagonal approximation requires that 

^ji + = ^ja + ^i4) this means that Ij-^ = Ij^ and Ij,^ = Ij^ or that Ij^ = Ij^ and Ij^ = Ij^. Both cases 

result in identical contributions (the fraction of orbits where = is negligible). Using (11) and 

replacing the sums by integrals with the density of diffracting segments jj resulting of multiplication 

of ( p!5D by 4, to take account of the various possible parities, leads to 

One should note that there is a non diagonal contribution from combinations of two twice diffracting 
orbits but this contribution is of higher order in r. 

The contribution which results from combinations of once and twice diffracting orbits is computed 
along the same lines as that from the combinations from twice diffracting and periodic orbits which 
was computed earlier. Using (|TD|) and (|TT|), the relevant contribution to the correlation function is 
found to be given by 

R^^^''\r]) = (j2 (-1 )^.i+^-^.i+^.2+^^.2+^.3+^^.3 V^^^^'^ ^ ^^'^^ n/}|2/)*e^'=('n-'.2-'.3) 



X e-*^('^i+'^2+'.3)e^¥ + c.c. + -v)] ) A' (3^ 
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where ji denotes the once diffracting orbit with length Ij-^ and j2,j3 denotes the twice diffracting 
orbit with length Ij,^ + /j3- Only combinations of orbits for which Ij-^ ~ Ij^ survive the smoothing. 
This leads to a saddle manifold of orbits that satisfy Nj-^ = Nj^ + Nj^ and Mj-^ = Mj^ + Mj^ (note 
that on this saddle manifold Nj-^ + Nj^ + Nj.^ + Mj^ + Mj^ + Mj^ is even) . To compute this contribution 
the sum over ^25^3 is restricted to the saddle manifold, on which the length difference is 

hi ~ ^2 ~ hi — ~2/^ ~ 2/3 ^^^'^ 

where ^2 + ^3 = ^ji, in analogy with (|3TD . The summation on the saddle manifold is then replaced 
by integration over t, I2 with density 



DpD*e"'K^+^)e-w'«+*¥ +C.C.+ iv^-v) 



3 



ll 



The contribution to the form factor is obtained by Fourier transforming (^D|) and replacing the sum 
over ji by an integral with the density resulting in 

The last contribution to be computed is from combinations of periodic orbits and orbits with three 
diffractions. The calculation is very similar to the one just performed. Using (^) and ([T2|) the 
contribution to the correlation function from combinations of periodic orbits and orbits with three 
diffractions is found to be 



(42) 



where the factor (— l)^Ji+^i2+^J3+*-^n+*^^2+^i3 was omitted since on the saddle manifold it reduces 
to unity. The saddle manifold is composed of orbits which satisfy 

Mj, + Mj, + Mj, = Mp. (43) 

The dominant contributions to the correlation function are from diffracting segments which are 
almost parallel to the periodic orbit as is shown in Fig. H. 
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FIG. 4. The configuration of lengths on the saddle manifold composed of periodic orbits and orbits with three difi^ractions 



The length difference is given by 



Ip Ij^ 



I- ~ 

n — 



Itj \tl 



2 



(44) 



2, li 2/3 2 I2 

The summation over ji,j2;j3 in (F^I) is restricted to the saddle manifold and then replaced by 
integration over /i, ti, Z2, ^2 where li varies between and Ip and I2 between and Ip — li. When the 
sums over ji, ^2 are replaced by integrals a factor of resulting from the density of contributions 
is introduced. The integral over ti,t2 is Gaussian and can be evaluated 



dtidt2 e 



*2 I (tl-t2) 



27r 

ik \^ /. 



p 



The resulting contribution to the correlation function is given by 



E 



dli 

JO 



dli 



-i{D^)*e~'^^'' + c.c + (77 -> -77) 



(45) 



(46) 



The integrals can be computed and then 
bution to the form factor 



is Fourier transformed leading to the following contri- 

(47) 



A 



■ 



The sum over periodic orbits can be replaced by an integral using the density ([T5|) , leading to 



(48) 
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Combining all the different contributions, that are given by (pj]), (|35D, (|37D, (^T]) and (|48|), the 
form factor is 

K(r) = 1 + '^Dt + - 3?/}^) + + 2|D|253Dr=^ - ^^^DV + O(r^). (49) 

The diffraction is assumed to conserve probability and therefore the diffraction constant has to 
satisfy the optical theorem. For an angle independent scatterer in two dimensions it takes the form: 



\D 



2 



= (50) 
leading to 

|D|4 = 8(|D|'-3?D'). (51) 



With the help of the last equations (|49| ) can be simplified to 



K{t) = l-^r + V + yD\* - -\D\'j + O(r^). (52) 

The term linear in time results from the combination of (forward) diffracting orbits and periodic 
orbits and is related, by the optical theorem, to the total cross section of the scatterer and thus 
its sign is always negative. The resulting form factor is 1 at r = and decreases with r due to 
the scatterer. At larger r the form factor starts to rise, and will get back to 1 for large r since 
the spectrum is discrete . Qualitatively, this type of behavior should be observed for integrable 



systems with point perturbations, since the forward diffracting orbits are always on a family of 
periodic orbits. From the terms to order one may get the wrong impression that the form factor 
depends on D and r only via the combination \D\'^t and in particular its value at the minimum is 
independent of \D\. This is not correct as one finds from the term of order r^. In particular the 
expansion exhibits dependence on |Dp of the minimum that is found reasonable compared to the 
numerical calculations presented in Sec. |V[ For this reason the calculation was terminated at the 
order r^. We now turn to the case where the scatterer is located at some typical (irrational) position. 
Only slight modifications of the calculation are needed, and these are related to the lifting of length 
degeneracies. 

B. Scatterer at a typical position 

When the scatterer is in the center of the rectangle the lengths of all the diffracting orbits with 
the same indices Ni were the same. If the scatterer is not in the center this is no longer true. An 
example is given in Fig. ^ 
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FIG. 5. All (1, 1) diffracting orbits for various scatterer positions 



The four dashed hnes in Fig. |^ represent the four (1, 1) once diffracting orbits which go from 
the scatterer to a corner and then back, they all have equal lengths. When the scatterer is not 
in the center the orbits with the same index {Nj,Mj) (solid lines) have different lengths. This 
situation is typical for odd-odd orbits (orbits where both Nj and Mj are odd). If the scatterer is 
shifted by {ttxSx, CLy6y) from the center the lengths of the odd-odd (A^^-, Mj) orbits take all four values 
= ^al{Nj ± 26xy + aliMj ± 25^)^. It will be assumed in what follows that 6x and 6y are typical 
irrational numbers, therefore the various lengths are not related in any simple way. For even-odd 



orbits the lengths take the two values = Ja^Nf + aii^Mj ± 25yY while for odd-even orbits the 



lengths take the two values Z^"''^'' = ^al{Nj ± 25^.)^ + a^M|. From (|A17|) one sees that these are 
related by time reversal symmetry. For even-even orbits all lengths take the same value /y^'*^^ = 



Nfal + Nja'^ that is identical to the value found for the periodic orbit {Np, Mp) = {Nj/2, Mj/2). 
Examples of such (unfolded) orbits are presented in Fig. |. 
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FIG. 6. Once diffracting orbits for a typical scatterer location. The solid line represents the (1, 3) odd-odd orbit, the dashed 
lines represent two of the even-odd (2,3) orbits, the dashed-dotted lines represent the odd-even (1,4) orbits and the dotted 
lines represents the even-even (2, 4) orbits. 



The four orbits in Fig. |^ that are in the first quadrant have all /i = 1, this notation is introduced 
since orbits in different quadrants generally have different lengths. Alternatively one can denote 
orbits using negative values of their index. The value of /i (see ( |A18|) and ( |A25|) ) thus determines 
the signs (of 5^ and 5y) in the formulas for the lengths of the orbits. The dashed line orbit in the 
second quadrant in Fig. || is also a (2, 3) orbits but with fi = 2. Its length is equal to the length of 
the (2, 3) orbit with /i = 1. This results from the fact these orbits are related by time reversal. The 
other (2, 3) orbits (/i = 3, 4) which are not shown have the same length which differs from that of the 
fi = 1,2 orbits. This situation is typical to all even-odd and odd-even orbits (the latter have fi = 1,4 
and are marked by dotted dashed lines in Fig. ^). The lengths of odd-odd orbits are different for 
different fi (these orbits are their own time reversals). The lengths of even-even orbits, marked by 
the dotted lines in Fig. ^, do not depend on fi and are the same as the lengths of periodic orbits. 

For this reason the contribution of the even-even diffracting orbits with a single diffraction com- 
bined with one of the periodic orbits is identical to the one obtained when the scatterer is in the center 
and reduces to K'^^'^\t) of (|16|) . Also the other contributions to the form factor are simply related to 
the ones found if the scatterer is at the center. For odd-even and even-odd {Nj, Mj) the lengths are 
divided into two degenerate lengths and therefore their amplitude is A^^^ /2 (the sums over /i in (|A25|) 
consists of two pairs of identical terms) and their density is 2p{lj) (see (p!0|) and ([TsD). The reason is 
that the number of orbits of each length is reduced by a factor of 2 and the number of values of lengths 
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in some interval increases by a factor of 2. Since the form factor is proportional to \ Aj\ and to p(lj) 
the resulting contributions to the form factor are ^K^^'°~^\t) = ^K^^''^~°\t) = ^'^g'^ (see (p!?!)). For 
similar considerations (all terms in the sum over fi in ( |A25|) are different) for the odd-odd orbits the 
amplitude is \A^^^ and the density is 4p(/j) leading to the contribution ^K''^'°^°\t) = ^|Dpr^ to 
the form factor. The contribution to the form factor corresponding to (|T8|) is 

ir(M(^) = ^\D\^r\ (53) 
16 

Therefore if the scatterer is located at a typical point the total contribution to the form factor, 
resulting from periodic orbits and once diffracting orbits corresponding to (|19|) is 

ir«(r) = l + !3Dr+ — |Z}|V. (54) 

16 

We turn now to look for the saddle manifold when the scatterer is at a typical position {axS^, CLySy). 
In order to find a saddle manifold one expands the length of the diffracting segments as was done 
in (|2T|). Since the shifts are of order of unity they are combined with 6N or 6M by defining 

5Nj =SNj±5x 

5Mj = 5Mj ± 6y. (55) 

The appropriate sign is the sign that appears in the length that is expanded. Using these definitions 
the rest of the calculation is identical to the one performed for the scatterer at the center except 
that the expansion is done with respect to 6N, 6M. The condition ( PTj) is replaced by 

= (56) 
Since {6x, Sy) are irrational this condition can be satisfied only if the signs of the shifts of both orbits 



are opposite and also if the condition (27) is satisfied. Note that the conditions for two diffracting 
orbits to be on the saddle manifold force these orbits to have the same parity. The contribution of 
each type of orbits to the form factor is calculated in what follows. For orbits of the odd-odd type 
both the shifts appear in the length, thus if the orbit ji has /i = 1 (the shifts have negative signs) 
then to satisfy (^) the j2 orbit has to have z/ = 3. Similarly, if one orbit has /i = 2 then the other 
one has to be of the z/ = 4 type. Therefore, from the 16 combinations that appear in ( |A25| ) only 
four will satisfy (pB]). The density of these contribution is ^ since the odd-odd orbits are a quarter 



of all the orbits. The calculation of the contribution is exactly as was done in ( p2[]3^ ), however the 



two factors of | result in 
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K'(^'°-"\t) = -—T^^D\ (57) 
16 

The contribution from the odd-even and even-odd orbits is similar, the only difference results from 
the fact that one of the shifts does not appear in the length of the orbit. This results in a degeneracy 
of lengths for orbits with different fi. Thus if one orbit has /i = 1 then not only the orbits with u = 3 
will contribute but also the orbit with u = 2 (or z/ = 4) marked by the dashed line (dashed-dotted 
line) in Fig. § which has identical length. Therefore, there are 8 combinations of n, v that will 
contribute to the saddle manifold leading to 



The contribution of the even-even orbits is 

K''^^'^-'\t) = --^t^^D\ (59) 

which is exactly as it was for a scatterer at the center. This is a result of the fact that the length of 
these orbits does not change when the scatterer position is changed. The sum of the non diagonal 
contributions of order is 

K'^^){T) = —^{-D^y. (60) 
16 

The computation of the terms is fairly similar to the case in which the scatterer is at the center of 
the rectangle. Therefore, only the relative factors resulting from the lifting of degeneracies of lengths 
will be calculated in what follows. The twice diffracting orbits are composed of two segments ji, 
j2- In the diagonal approximation each segment is summed independently. For the diagonal terms 
resulting from once diffracting orbits it was shown that a relative factor of ^ is obtained in ( ^41) 
compared to ([T9|). Since the sums over ji and j2 are independent a relative factor of (^-^ = ^ 
is obtained. The diagonal contribution of twice diffracting orbits for a scatterer at some typical 
location, corresponding to (^Tl), is therefore, 

K^'\r) = ^^\D\^r\ (61) 

We turn to compute the relative factor that results from relocating the scatterer from the center to 
a typical location for non diagonal terms of once and twice diffracting orbits. For this contribution, 
the saddle condition ( p6D is replaced by 

6Mj, =6Mj,+6Mj,. (62) 
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To understand the effect of the location of the scatterer the non integer part of (|62|) is studied. 
This non integer part depends on the displacements of the scatterer from the center Sy. It is 
convenient to discuss different parities of ji separately, giving a factor of | to each case (since each 
case contributes a ^ of the total contribution when the scatterer is at the center). If the orbit ji 
is even-even then the parity of j2 should be equal to the one of js. This is exactly identical to the 
calculation of the non diagonal contribution of order which resulted in a relative factor of ^ 
in ( [5D| ) compared to (^). Therefore, the contribution when ji is even-even, is ^ of the total term 
for the scatterer at the center, given by (^l]) . 

If the orbit ji is of the even-odd parity then the combinations of j2 and js that satisfy ( |62|) are: (a) 
j2 even-even; js even-odd and (b) j2 odd-odd; ja odd-even and also two combinations where the roles 
of j2 and j3 are interchanged. In case (a) there are 8 combinations of fij^, fij^ that will satisfy ( |62D 
for each /ij^ (compared to 16 when the scatterer is at the center). Since this case (with ji even-odd) 
accounts for a | of all contributions it results in a relative factor of = |. In case (b) there are 4 
combinations of /i^j and ^j.^ that satisfy (|62D for each /i^^. This results in a relative factor of ^. The 
relative contribution of terms where ji is even-odd is thus given by | (^2| + = ^, where the 
factors of 2 result from the cases where the parities of j2 and js are interchanged. The case where 
ji is odd-even gives an identical contribution to the case in which ji is even-odd. 

When ji is odd-odd (^) is solved if (c) j2 is even-even and js is odd-odd or if (d) j2 is even-odd and 
j3 is odd-even. There are also two additional solutions that result from interchanging the parities of 
j2 and js- One can verify that for all the parities of j2 there are 4 such solutions for each value of 
/iji compared to 16 solutions when the scatterer is at the center. This results in a relative factor of 
ll = Yq compared to (|l]). 

Combining all the contributions from different parities of ji one finds that for a scatterer at a 
typical location the contribution of non diagonal terms from once and twice diffracting orbits is 
^ + 2^+j|^ = II relative to the case where the scatterer is at the center. Therefore using (0) one 
finds 

i^'(3a)(^) ^ '^\D\^qDt^ (63) 

The contributions from non diagonal terms of periodic orbits and orbits with three diffractions 
should satisfy the condition 

E SN,, = 

i=l 

E SM,^ = (64) 

i=l 
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corresponding to (|56D. The non integer part of this condition behaves exactly in the same way as 
the non integer part of (|6^). This is true since both signs of the displacements always appear. This 
results in the fact that the factors relative to the case where the scatterer is at the center are equal 
for both of these terms. Therefore using ( ^S]) one finds 

25. 



T 



96 



(65) 



Adding all the contributions 
position of the scatterer 

9 



pip, ( p5D and ( p^D leads to the form factor for a typical 
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K{t) = 1 + '^Dt + i-\D\^T^ + —^{-D^y 
16 16 



512' 



4„3 



32' ' 96 ^ ^ 



The optical theorem 

K{t) = 1 



and (Elh can be used to obtain 



-r + 
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\D\'t' 
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D|V + 0(r^). 



4 128' ' 512' ' 1536' 
This form factor has the same qualitative features of the form factor ([5^), it is 1 for r 
decreases linearly for small r. See the comment following (0). 



(66) 

(67) 
and 



III. A MODEL OF AN INTEGRABLE BILLIARD WITH A POINT SCATTERER 



An easily solvable model is that of a point interaction [^,^. This interaction is the self-adjoint 
extension of a Hamiltonian where a point is removed from the domain. It can be viewed as formally 
represented by a 5-function potential. The influence of such an interaction on the spectral statistics 
was first investigated by Seba [|lT|. The system that was investigated was the rectangular billiard 
with a point interaction, and the spectral statistics were found to differ from Poisson statistics when 
the point scatterer is added. The spectral statistics of this system and of similar systems were 
the subject of many works [r^j36|j3^j4^j4^j4^j45|j46|j4^j48|] . Since the energy levels of this system are 
obtained from the zeroes of a function (as will be explained in the following) it is possible to compute 
a large number of levels and thus also the form factor with considerable accuracy. Therefore, we use 
the rectangular billiard with a point scatterer to obtain numerically the form factor and compare it 
to the predictions of the analytical theory, presented in the previous section. 

Point interactions are useful since the Green's function of the problem G^{z; x, x') can be expressed 
by the Green's functions of the problem in absence of the scatterer G'(z;x, x') 



G^{z; X, x') = G{z; x, x') + T{z, ^^'(-z; x, :>Co)G{z; Xq, x') 



(6^ 
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where ^ is the self adjoint parameter that is related to the interaction strength, xq is the scatterer 
location and Q^/z > 0. T is the on shell transition matrix and it is given by 

-1 



T{z,0 = {e^^-l) |dyG'(z; y, xo)G'(M y, xo) + e*«(Az + z) J dyG{z;y,^o)Gi-At;y,^o) 



(69) 



A is an energy scale usually set to be unity, since there is only one parameter in the problem which is 
the interaction strength. The point scatterer is characterized by the self-adjoint parameter ^, while 
T is independent of angles, resulting in s-wave scattering. The energy levels of the perturbed system 
are the poles of the resolvent, and if they do not coincide with the eigenvalues of the unperturbed 
system they are also the poles of the transition matrix T. The unperturbed Green's function of the 
rectangular billiard is given by 

G(.;x,x') = £*i<=^ (70) 



0. (71) 



where ipn and En are the eigenfunctions and eigenvalues of the billiard. Substituting (^) into ( pQf ) 
leads to the equation for the energy levels of the perturbed problem 

In the geometrical theory of diffraction used in the previous section the effect of the scatterer that 
is located in the billiard was expressed by the diffraction constant D that determines the scattering 
in free space. The diffraction constant in free space is the T matrix ( |69D if G{z; x, x') is the Green's 
function in free space. It was calculated in App. B leading to 

D{z) = — , , . , . (72) 

2 "Ha; 4 l-cos^ ^ 2 

The self adjoint extension parameter in free space is and it need not to be the same as ^ that was 
used for the bound problem (|7lD. The reason is that the self adjoint extensions depend not only 
on the scatterer but also on the boundary conditions. The general relation between extensions of 
different problems is unexplored to the best of our knowledge, however, a connection can be made 
in the case that is considered here. 

The energy scale associated with the scatterer is A. For A ^ A the spectrum is effectively contin- 
uous, therefore the T matrix of the bound problem is approximately the diffraction constant of the 
unbound problem, and also the self adjoint extension parameters ^ and should be approximately 
equal. Another way to establish this correspondence is to compare length scales. The length scale 
associated with the scatterer is l/y/A and the condition A A is that l/VX <^ \CA. where A 
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is the area of the bilhard. The condition is that the length associated with the scatterer is much 
smaller then the dimensions of the billiard. This is also the condition for the validity of the semi- 
classical geometrical theory of diffraction used in the previous section. The scattering is strongest 
{\D\^ maximal) when D = —4i. The variation of D{z) with z is weak and therefore its value at the 
maximum is expected to be a good approximation for a whole energy region around the maximum. 
Substitution of D = —Ai in ( |^ ) and ( |B7D yields: 

1 9S 

K{t) = 1-At + ?>2t^-—t^ + 0{t^) (73) 

3 

when the scatterer is at the center of the rectangle and 

1 ^7 

K{t) = 1-At+1^t^ -=^T^ + 0{t^) (74) 

6 

when the scatterer is at some typical location. 

When the scatterer is at the center of the rectangle the system has a symmetry. The (unperturbed) 
eigenfunctions that vanish on the scatterer remain eigenfunctions of the perturbed problem, and 
so are their eigenvalues. To obtain the other eigenvalues one has to use only the non-vanishing 
eigenfunctions in (ffll). These are the eigenfunctions that are even functions with respect to reflection 
with respect to both the X and Y axes that pass through the center of the rectangle. Since the value 
of these eigenfunctions for all eigenvalues is the same ('?/'n(x = 0) = -^=) at the center of the rectangle, 
the resulting equation (^) is the same as for the Seba billiard with periodic boundary conditions ||13|| . 
The spectral statistics of the Seba billiard with periodic boundary conditions were also shown to 
be the same as the spectral statistics of some star graphs PB|JT^ . However, ( |75D describes the form 
factor of the full spectrum. Therefore, in addition to the eigenvalues that are perturbed due to the 
scatterer, all of the eigenvalues that are unaffected are included in the calculation of the form factor. 

The form factor of the perturbed levels is related to the form factor of the full spectrum. It is 
composed of the even-even eigenfunctions of the rectangle. The eigenvalues of the four different 
symmetry classes of the rectangle can be assumed to be uncorrelated, therefore their combined 
form factor is the sum of the form factors (weighted using their density) of the different symmetry 
classes, with r rescaled as required by (|I|) and (H). Since the unperturbed levels have Poisson level 
statistics and since levels belonging to the four symmetry classes have the same mean density (in 
the semiclassical limit) the relation between the full form factor and the form factor obtained from 
only perturbed levels is given by 

Kf.,u{T) = \ + \K,,r{^T). (75) 

The perturbed form factor is exactly the form factor that was computed for the Seba billiard with 
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periodic boundary conditions and for star graphs. The first powers in r of the form factor for the 
star graphs are given by 

KpUr) = 1 - 4r + 8r2 - ^r' + ©(r^), (76) 



as can be seen, for example, from equation (14) of ||3^ where the form factor is given up to order of 
r*. Using (|75D gives exactly the form factor (|73|). 

IV. NUMERICAL RESULTS 

The spectral form factor is not easy to compute numerically. The form factor is a very oscillatory 
function of time and much averaging is needed in order to obtain its mean behavior. This requires 
computation of many levels, which is not always possible numerically. However, such calculations 
are relatively easy for the point scatterer defined in the previous section by equations (^) and (|69|). 
The eigenvalues can be calculated from ([7T|) . This equation in particularly convenient since J2n z-e 
is monotonically decreasing when z is varied between any adjacent unperturbed levels, ensuring that 
every level must be in an interval between two adjacent unperturbed levels En and -En+i (unless 
the unperturbed wavefunctions vanish at the scatterer). To compute the energy levels numerically 
equation ([7T| ) is solved, for example by truncating the sums at some large n and replacing the tails 
by integrals. 

This equation is solved in some energy window. The energy scale parameter is chosen so that 
A ^ A, and therefore D and ^ can be approximated by their values for a scatterer in free space, 
as discussed in the previous section. In calculation presented here ^ was chosen so that D = — 4i, 
resulting in maximal scattering, for z in the center of the window. The other values of D that were 
also examined, are D = 2 — 2i and D = |(2 — i). Since D does not change much with z the results can 
be compared with the analytical predictions obtained within the semiclassical geometrical theory of 
diffraction. 

For a scatterer at a typical location the form factor is presented in Fig. ^. 



25 




0.0 0.5 1.0 1.5 2.0 

X 

FIG. 7. The form factor of a scatterer at a typical position, for some diffraction constants, thin lines, compared to the 
analytical result (^^, heavy lines. 

The form factor is computed from energy levels between 20000 and 40000. The mean level spacing 
is A = 1 while A = 1000. To smooth the form factor an ensemble average over 20 systems was 
performed. The aspect-ratio of each systems was taken as a random number between 1/2 and 
1. Additional smoothing was performed by averaging the results found for 4 consecutive values of 
r. The position of the scatterer was also chosen at random, but it was the same position for all 
ensembles of aspect ratios. Similar results were found for several other scatterer positions. The 
agreement between the numerical results and the theoretical prediction ( |67D is good at short times, 
where the expansion truncated at a low order is a good approximation. 

For the scatterer at the center (xq = 0) only levels with eigenfunctions that are symmetric with 
respect to the X and Y axes are perturbed by the scatterer as is clear from (^). Based on the 
assumption that the spectra of the various symmetry classes of levels are uncorrelated (|75D was 
derived. To test it the form factor of the full spectrum and of the levels which are affected by the 
scatterer were calculated. These are shown in Fig. |. The calculation is done for parameters similar 
to the ones used in the case where the scatterer is not at the center, and an averaging over 70 systems 
is performed here, while the energy range considered was 60000 — 80000. 
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FIG. 8. The form factor of a scatterer at the center (xo = 0), for only perturbed levels as well as the full form factor 

The form factor of the rectangular bilhard was included for comparison, demonstrating that at 
short times all form factors deviate from 1. This deviation is due to the fact that one examines times 
of the order of the period of the shortest orbits and there are not enough orbits that contribute to 
the form factor for such times. In order to demonstrate it, the form factor of the rectangle was 
calculated from the exact energy levels in the same way as the form factors with diffraction. The 
result is shown in Fig. || and this form factor is also found to deviate from 1 at short times. The 
form factor calculated from all levels and the scaled form factor obtained from the perturbed levels 
with the help of ( |75D are compared to the analytical result ([75| ) in Fig. ^. 
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FIG. 9. The form factor for the scatterer at the center (xo — 0) compared to the scaled form factor calculated from perturbed 
levels by ( [ts] ) and the analytical result ( |7^ ) 



It is clear that the symmetry argument (|75D is vahd. The scaled form factor is indeed the same 
as the one that is obtained from the full spectrum. Both slightly deviate from the analytical pre- 
diction ([T3|). The deviation at short times of the full form factor is due to the fact that there are 
not enough relevant short orbits as was discussed earlier, and at these times the scaled form factor 
is probably more accurate. Since the periodic orbits prediction is in agreement with the form factor 
studied in [P!^PU| , P7| (that is exact) the deviation probably results from the numerical solution, or 
of the fact that the energy is not high enough. 

The agreement between the analytical results obtained in Sec. Jl| and the numerical results of the 



model of the point scatterer outlined in Sec. |T| was found for A ^ A, where for the point scatterer 
the values of ^ and D of a scatterer in free space were used. The requirement A ^ A is essential 
for this agreement. For A ^ A if one uses D = —4i and the corresponding values of the self adjoint 
extension parameter ^, found for the scatterer in free space, the effect of the point scatterer on the 
spectrum of the bounded system is even not maximal (as is the case when A ^ A) and clearly there 
is no agreement with predictions of Sec. for these values of A and D. 



V. SCATTERING FROM A LOCALIZED POTENTIAL 



So far it was assumed that the diffraction constant is given. In this section the diffraction constant 
for a given potential, which is localized around some point, is calculated. Consider a potential which 
is assumed to have a length scale a such that 
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U{r) = (-] (77) 



where /(y) is small for y that is large compared to unity. The scattering from such a potential at 
energies for which the wave number k satisfies, fca -C 1, is studied in what follows. 
The diffraction constant is the on shell matrix element of the T matrix 

D{e\e) = {k\T{E)\q) (78) 

where k is the outgoing momentum (in direction 6') and q is the incoming momentum (in the 
direction 6). The energies of the incoming and outgoing waves are equal, that is 

k = q = ^ (79) 

(in units h = 1 and m = ^ used in the paper). In order to compute the diffraction constant it is 
convenient to start from the Lippman-Schwinger equation for the T matrix 

T{E) =gV + gVG{E)T{E) (80) 

where G{E) is the Green's function for free propagation at energy E and g is an expansion parameter 
which will be set to unity at the end of the calculation. One might be tempted to compute the on 
shell matrix elements of T using the Born approximation which is obtained by iterating equation (ROl) 



(k|T(^)|q) = ^(k|U|q) + /(k|UGU|q) + ^3(k|UGUGU|q) + ■ • • (81) 

However, in the limit ka —>■ the Green's function diverges as In ka and therefore the Born series 
diverges. 

A method which behaves regularly at low energies {ka ^ 1) was developed by Noyce |]50|. In this 
method the scattering in the forward direction (another chosen direction can also be used) is in some 
sense resummed. In two dimensions this method was used to describe low energy bound states 



and low energy scattering [^. The Noyce method is derived in the following. The T matrix is 
written as a product of the diagonal matrix element times a function of the incoming and outgoing 
momenta, 

(k|T|q) = Mk,q)(q|T|q). (82) 
Taking the diagonal matrix elements of the Lippmann-Schwinger equation (RDf) leads to 



/ IT-I \(^ ^ f ^^1 (q|U|ki)/i(ki,q)\ 

where the Green's function in the momentum representation is 
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(27r)25 (ki -k; 
E-kl + ie 

Note that ki is not on the energy shelL 

Using (|83|) in the non diagonal matrix elements of ( ^O]) leads to an integral equation for h: 

(ki|U|q)(q|U|k2)\ 



(q|U|q) (27r)^ J \ 



(q|U|q) 



^ /i(k2,q). (85) 



+ ie - 

This equation can be iterated resulting in an expansion of h in powers of g. Therefore with the help 
of ( p2D and ( ]83| ) the T matrix can be written in the form 

(k|T|q) = . (86) 

Both /i(k, q) and the denominator -Dg(q) can be expanded in powers of g as 

(k|U|q) / (k|UGU|q) (k|U|q)(q|UGU|q) \ 
^^^'^^ (q|U|q)+H (q|U|q) (q|U|q)2 j 

"(k|UGUGU|q) (k|UGU|q)(q|UGU|q) 



(q|U|q) (q|U|q)2 
(k|U|q)(q|UGUGU|q) (k|U|q) (q|UGU|q)2^ 



(q|U|q)2 (q|U|q)3 



+ ••• (87) 



and 



(q|UGU|q) 2^(q|UGUGU|q) (q|UGU|q)^ 

^<;(q) = 1 / iTTi \ 9 



(q|U|q) " V (q|U|q) (q|U|q)2 J 

■(q|UGUGUGU|q) ^(q|UGU|q)(q|UGUGU|q) 



(q|U|q) (q|U|q)2 

,88) 

While the expansions for h(k,q) and -Dg(q) seem to be complicated they are actually composed 
of all the ways to break the matrix elements, containing chains of the form UGUG ■ ■ ■ UGU with 
several Green's functions, to products of smaller elements with less Green functions. Each order of 
g consists of terms with the same number of Green's functions in it, and when an element is broken 
into a product an extra factor of ^q|"^^|q^ is added. In the expansion for -Dg(q) all factors contain 
at least one Green's function but in the expansion for /;,(k, q) factors of the type (k|U|q) appear. 
Consequently the g^ term for /i(k, q) can be explicitly written. It includes 2^ terms (and 2^~^ 
terms for Dg{q), some of which may be identical). Note that jy^j^ can be expanded in powers of g 
leading back to the Born series (|8lD. 
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To compute the diffraction constant the matrix elements are computed and the length scale a of 
the potential is scaled out. For example 



(k|UGU|q) 



d''ndW''^''Uir,)Goir^,r2)Uiv2)e'^■'^ 



^9) 



where Yi = ^ and ( [77D was used and Go is the free Green's function in the position representation, 
namely Go(ri,r2) = (ri|G(-E)|r2) = Go(A;a; yi, y2), where k = ^JE. These matrix elements are 
expanded for ka keeping terms of order k^a^ . Then the series in the parameter g is summed. 
Both the exponentials (k|r) = e"*'' '" and the Green's function should be expanded to the same 



order. In the limit of small ka one finds 

Go(fca;yi,y2) = 

where 
and 
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0(yi,y2) (90) 



(91) 



0(yi,y2) = TT^^I^i ~ ^2! 

ZTT 



(92) 



while 7 is Euler's constant. The leading order part of the Green's function which is coordinate 
independent is d\ and, as was mentioned earlier, it diverges logarithmically when ka 0. In App. C 
it is shown that in many of the terms to be computed this constant part will be canceled. This is 
the advantage of the Noyce method compared to the Born expansion for the present problem. 

Calculations presented in App. C lead to the asymptotic expression for the angular dependence of 
the diffraction constant 



D( 



V{g)-gVod, + k^a^Q{g) 
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where c.c. stands for complex conjugate. 



k=l 



Mn 



V\g)Hg)\' 



(93) 



(94) 



(95) 



31 



and 



M++ = {V{g) - gd.Vo) ^— + V{g)gd,Vo^ (96) 
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M+_ = - gd.Vo) ^^^^ + (97) 



W{g)\ 
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M_+ = {V{g) - gd.Vo) ^^—^ + V{g)gdiVo^^ (98) 



M__ = {V{g) - gd.Vo) + V{g)gd,Vo'-^^i^. (99) 



■ ■ ----- \/o2 

The constants Vo,V{g),n^"''''\g),a{g),a{g),S{g), and £fc are angle and energy independent and are 
given by a series of integrals over the potential. They are defined in (P2|) , (P^D , and ( |(J14| )-( [C2(]| ). 

Equation ( PBj ) is the main result of this section. It provides the functional form of the diffraction 
constant for ka <^ 1. It also gives a recipe for the computation of the coefficients of the angle 
dependent factors as sums of certain integrals involving the potential (^). The physical scattering 
constant is obtained when g = 1 is substituted in (^31) . It is of interest to note that the simple 
dependence on the angles 6 and 6' is expected since formally it should be the same dependence as 
in the Born series ( |8TD which is also simple. This diffraction constant will be used in Sec. ^for the 
calculation of the effects of angular dependence on the form factor. 

Finally, one can verify that the diffraction constant (pSf) satisfies the optical theorem 



QD{9,9) = --^ r^rf^' \D{9',9)\^ (100) 



TC Jo 

2^2 



up to order k a . 

VI. THE EFFECT OF ANGULAR DEPENDENCE ON THE FORM FACTOR 

In this section the form factor resulting of angle dependent diffraction is calculated to order r^. 
The complicated parts of the calculation are not related to the dependence on the angles but rather 
to the contribution from non diagonal terms that were calculated in Sec. |l|, therefore this calculation 
will not be repeated here, and only the modifications needed to describe angle dependent diffraction 
constants will be presented. 

Since the diffraction constant depends on angles, the density of diffracting orbits of length / and 
outgoing direction 9 is needed (the incoming direction depends on the outgoing direction, via (|A17| )). 



For simplicity consider first the case where the scatterer is at the center. Instead of orbits that leave 
the scatterer, hit the walls, and return to the scatterer it is possible to examine orbits which connect 
the scatterer to one of its images, when the plane is tiled with rectangles, as described in detail 
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in App. A. These orbits have the same length and the same outgoing directions as the diffracting 
orbits. To compute the (averaged) density of such orbits the number of orbits with lengths in the 



range (1 — ^,1 + ^) and outgoing directions in the range (6^ — y, 6* + y) is estimated. The area of this 
domain is 16661. When computing the form factor, in the semiclassical limit, the lengths of orbits 
that contribute scales as k. There is a natural width 61 which results from the averaging in (0) and 
is taken such that 1^5/^/. Therefore, for any small but finite 66 the area 16661 is large in the 
semiclassical limit. The number of orbits (of a each type) in such area is given by jj666l. If the 
diffraction constant depends slowly on the angles, that is, its value changes considerably only when 
the angles change by 56* ~ 1 (namely 0{k^)) then the sums over diffracting orbits can be replaced 
by integrals with the density 

pa^) = ^- (101) 

On this scale the angles of diffracting orbits are distributed uniformly. Note that if the diffraction 
constant is angle independent one can integrate ( |101|) over 6 in the interval (0, |) and obtain (p!5|). 



As a result the form factor for an angle dependent diffraction constant can be computed if instead 
of the constant D of Sec. |I|one uses D{6' , 6) of (|93D with g = I with appropriate averages on angles. 
These averages depend on the type of the diffracting orbits that contribute and are computed in the 
following. 

A. Scatterer at the center 

To compute the contributions to the form factor each type of orbits should be considered separately. 
This results from the fact that the incoming direction is related to the outgoing direction in a different 
way characterized by the parity and by the indices fi and u of App. A, where it is explicitly given 
by ( A17 ). First the diagonal contribution from diffracting orbits is computed and later the non 



diagonal contributions from twice diffracting orbits and periodic orbits will be computed. 

First consider the contribution of the even-even orbits following the lines of the calculation leading 
to (|16|) of Sec. |T|. This is the contribution of forward diffracting orbits. For these orbits the incoming 
direction is the same as the outgoing direction. The index fi of ( |A25| ) specifies the quadrant of the 
angle. In Sec. H the factor resulting of summation over fi was included in the amplitude but for the 
angle dependent D{6',6) the averaging over angles is performed with the help of the density ( p.01| ). 
The summation over /i results in extension of the integration interval from (0, |) to (0,27r). When 
the scatterer is at the center orbits with different fi (but with the same index Ni, N2) have identical 
lengths. Since 
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and 
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the contribution of even-even orbits, that are the forward diffracting orbits, to the form factor is 
given by 
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This contribution corresponds to (p!6| ) of Sec. The diagonal contributions from the other types of 
diffracting orbits are computed in a similar manner and presented in App. p 1| . Substitution of the 
diffraction constant (031) in the diagonal contribution (p4|) leads to 



K«(r) = l-i|C|V+ ICpCV 



where 



C 



Vn 



V{1) - Vodi + Pa2g(l) 



and 



C' = l- ea^vHi 



ail) 



(105) 



(106) 



(107) 



See (0), (|94D, ([C^), (|C10| ), and ( |U20|) for the definitions of the various quantities used in (|106| ) 
and (11071) . 

An additional non diagonal contribution results from orbits which have almost identical lengths. 
At order such a contribution results from combinations of twice diffracting orbits and periodic 
orbits. For angle independent scattering this contribution was computed in Sec. y (see (|2^^3|) 
there). The modifications which result from the angle dependence are computed in App. p 2| leading 
to ( Pll|) that yields 



ir'(2)(r) = -(3?C2)C"r2 



(108) 
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where C and C are defined by ( |106| ) and (|107|) . To obtain the form factor of an angle dependent 



scatterer located at the center of a rectangular billiard, to order r^, equations ( |105|) and (|108|) are 
summed resulting in 

i^(r) = l-i|C|V + i|C|W. (109) 

This form factor resembles the form factor (^2]) that was obtained for angle independent scattering. 
It starts at 1 and decreases for small times until at some point it has a minimum. At larger times it 
rises and approaches 1 for long times. The angle dependence just modifies the constants in front the 
powers of r, but these changes are of the order k'^a? and typically cannot change the sign of these 
coefficients. 

B. Scatterer at a typical position 

When the scatterer is not at the center the lengths of the diffracting segments of the same parity 
vary and some length degeneracies are lifted, as discussed at the end of Sec. ^ The form factor for 
an angle dependent scatterer in some typical position is calculated along the lines of the calculation 
of Sec. |T[ Only the modifications due to the angle dependence are presented. First the diagonal 
contribution from once diffracting orbits is computed and then the non diagonal contributions from 
combinations of twice diffracting orbits and periodic orbits are computed. 

To compute the diagonal contribution it is convenient to separate the orbits into types determined 
by the number of bounces from the walls (or the parity of their index Nj). For orbits of the even-even 
type the lengths of all 4 possible segments (with different /x) are identical and do not change when 
the scatterer changes position. Therefore their diagonal contribution is identical to the contribution 
that was obtained for a scatterer at the center and is given by (|104|) . The other contributions are 
calculated in App. p 3| leading to 

ir«(r) = 1 - ^|Cp + ^\CW^C[y (110) 

where 

"3ft/i(+-)(l) l^(l)|a(l)| 



C[ = ?, + k^a^ [V{1) - Vodi] 



;iii) 



and the various quantities are defined by (pi]), ( |106| ), (|C2|) , (|C10| ), and ( |C20|) . This is not the 



total contribution to the form factor (in order r^). There is a non diagonal contribution from twice 
diffracting orbits and periodic orbits which is computed in App. P 4| and is described in what follows. 



The non diagonal contributions are from twice diffracting orbits and periodic orbits with almost 
identical lengths. The condition for such combinations to contribute is given by (|56D. This condition 
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ensures that the two segments of the diffracting orbit are of the same type, as it was for the scatterer 
at the center. The fact that the lengths of the segments might depend on the index /i was discussed 
following the condition (^61). This discussion is also valid for angle depended scattering. As for the 
case where the scatterer is at the center the integral over angles can be factored out of the saddle 
manifold integral. 

For even-even orbits the lengths do not depend upon /i and the non diagonal contribution from 
these orbits and periodic orbits is given by ([IT^). The other contributions are calculated in App. |D 4| 
resulting in (|D19|) leading to 

K'^'\r) = -^MC'C[)r\ (112) 
Id 

where C[ is given by (|111|) . Finally summing (|11CI|) and (|112|) results in an expression for the form 
factor, correct to order r^, for the case where the scatterer is at a typical position, 

K{t) = 1~\\C\'t + ^^\C\'t\ (113) 



This form factor is rather similar to the one obtained for angle independent scattering (|67|). The angle 
dependence modifies the coefficients only slightly. The coefficient C satisfies the optical theorem (pO|), 
therefore this form factor is identical to one obtained for an angle independent potential, with 
diffraction constant C . 



VII. SUMMARY AND DISCUSSION 



In this paper the form factor K{t) was calculated for a rectangular billiard perturbed by a strong 
potential, that is confined to a region which is much smaller than the wavelength. In the first part 
(Sees, m, |ITl| and |rv| ) a highly idealized situation of a scatterer that is confined to a point and its 
action is represented by the diffraction constant that is a complex number like D = — 4i or 
D = 2 — 2i, was explored. The form factor K{t) was calculated to 0{t^) and the results are given 
by (0), if the scatterer is at the center, and by (|6^) if it is at a typical position. These results were 
compared to the ones obtained numerically for various values of D in Figs. |^ and ^. Reasonable 
agreement was found in the small r regime, up to the minimum. The numerical calculations were 
performed for the Seba billiard, that is a rectangular billiard with a 5-like scatterer inside. Its 
T matrix is ( |69D and the energy levels, used in the numerical calculation of the form factor, were 
obtained from the numerical solution of (^). This is an idealized model approximating a rectangular 
billiard with a strong potential that is confined to a very small region. The model is defined by the 
self-adjoint parameter ^ and the energy scale is set by A. The connection with the analytical results 
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of Sec. H is via the diffraction constant D of (|7^), that is defined in free space. If A ^ A, where A 
is the mean level spacing (or 1 / is much smaller then the dimensions of the billiard) , then the 
self-adjoint parameter ,^1 of ([T2|) is approximately equal to S,, that is the corresponding parameter 
in the presence of the boundary. This condition was found to be essential. The variation of D with 
energy is slow. For these reasons the numerical results of Figs. ^ and |^ are expected to agree with 
the analytical predictions of Sec. |I|, for small r, and indeed such agreement is found. 

In the second part of the paper (Sees. |V| and |VJ) a problem with an arbitrary perturbing potential 
U{v) of ([77|), that is large in a small region of extension a, was studied. For rectangles with such 
potentials the diffraction constant was calculated to the order (fca)^, where k = 2tx/\ = \/E is 
the wavenumber. The result is given by (0) with g = 1. The angle dependent terms are of order 
ka and {ka^. In the limit {ka^ ^ 1 or a <^ A the diffraction constant D{9',6) is approximately 
independent of angles and reduces to 

where Vq and V{g) are given by ( p2|) and ( |C9| ) while 7 is Euler's constant. One can easily see that 
Co satisfies the optical theorem. We do not know much about the convergence of the series (|C9| ) and 
( |C10|) for V{g). Therefore V{1) is assumed to be the analytic continuation from the small g region. 
Exploration of these series for various potentials is left for further research. If the potential is such 
that in the a — limit Cq is well defined, then in this limit it approaches the diffraction constant 
D that was studied in the first part of the paper, taking the constant value D in Sec. |I| and (ff^) in 
Sees, m and For this to hold it is required that in the limit a — also 

Vo 27r r V 2 / / 47r I A i 8 1 - cos^i 



If the function /(y) of ([77|) is independent of a, then also V{l)/Vo is independent of a, as is clear 
from ( |C2|) and (|C9D. The calculations of Sec. and App. C, leading to (0), do not require that / 
is independent of a. Assume for example that / factors in the form 

/(a,y) = /(^)/i(y) (116) 

where /i(y) is independent of a, while / is independent of y and /a is some characteristic length 
scale. If lima_»o /(q-Z^a) = then lima^o^(l) = 1- The factorization ( |116D implies 

^o = /(^) /A/i(y)- (117) 

Existence of the limit ( |115|) requires that 
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(118) 



where /o is a constant independent of a, in agreement with the scahng used in the mathematical 
hterature ( [0 p. 103). It leads to the natural identification l\ = j^. Therefore /a is the length 



scale associated with the energy scale A that was discussed after (72). The limit is characterized by 
the two parameters ^ and A. If the factorization ( |116| ) exists, then for a <^ A the results obtained in 
the first part of the paper provide a good approximation for the ones obtained in the second one, for 
arbitrary but well defined potentials U of (^). For any given potential f/(r) that is concentrated in 
a small region, of extension a, one can use (^) to compute the form factor. The condition for the 
applicability of the semiclassical approximation combined with the geometrical theory of diffraction 
(GTD) is 

a -C A <C a^;, ay, (119) 

where ax and the sides of the rectangle. 

If the scatterer is at a typical position then the form factor is given by ( |113|) . Note that C 
of Eq. ( |113| ) (defined by (|106|) ) satisfies the optical theorem, since QiX) is real. Therefore the 
form factor of the angle dependent scatterer reduces (to the order r^) to (|67D . That is, the angle 
dependent scattering is affecting the form factor in the same way (to the order computed here) as 
angle independent scattering. Therefore the angle dependence plays no role up to this order. For the 



scatterer at the center the situation is somewhat different. In this case, the form factor ( |109|) should 
be compared with (^). There is a correction C resulting of the angular dependence of D{6',6) 
given by (|9^) . This difference is a consequence of the increased number of length degeneracies of the 
diffracting orbits when the scatterer is at the center compared to the situation when it is located at 
a typical position. 

For the case where the scatterer is at a typical position, using angle dependent terms to order (ka)"^ 
does not change the result for the form factor, compared to the one found for the angle independent 
leading order. Therefore it is reasonable to assume that the results are robust and the limit a 
describes correctly the physics of the regime a < X. This is so although the classical dynamics (in 
the long time limit) are expected to be chaotic in nearly all of phase space and similar to the ones 
of the Sinai billiard. This improves the chances for the experimental realization of the results of the 
present work. Note that for a ^ A, semiclassical theory works and the system should behave as a 
Sinai billiard, with level statistics given by Random Matrix Theory (RMT), in some range [|7|,pO|. 

The spectral statistics found in the present work differ from the ones of the known universality 
classes. It is characterized by the form factor of the type presented in Figs. ^ and ^. A characteristic 
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feature of the form factor is that it is equal to 1 at r = 0, resulting of the fact that for small r 
the number of classical orbits that are scattered is small. The contribution that is first order in r 
originates from the term (14), that leads to the contribution QDt in (^^, (|49|) and (|66D. By the 



optical theorem ( |50D it is always negative. This is expected to hold for a larger class of systems, 
where there is forward diffraction along periodic orbits. For r ^ 1 the form factor approaches unity 
because of the discreteness of the spectrum |]3D|. These are the physical reasons for the qualitative 
form of K{t) depicted in Figs |^ and |^. 

The form factor was computed when the scatterer is at the center and when it is at a typical position 
shifted by {a^Sr^, aySy) from the center (with all numbers in this expression being irrational). The 
results are different since for the scatterer at the center there is a high degree of the degeneracy of 
the orbits involved. In this case the form factor is related by the symmetry argument (^) to the one 
found for periodic boundary conditions, that is known exactly. The validity of ( [75D is demonstrated 
in Fig. 1^. Off diagonal contributions of orbits belonging to saddle manifolds turn out to be of great 
importance (see (pOD and the calculations that follow). In view of the work of Bogomolny |^ this 
should be generic for integrable systems in presence of localized perturbations. 

A natural question that should be explored is whether the problem of the rectangular billiard with 
a point scatterer, that was studied here, represents a larger universality class and whether it can be 
related to some RMT models. 
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APPENDIX A: DIFFRACTING ORBITS CONTRIBUTIONS 

In this Appendix the contributions of diffracting orbits to the density of states are calculated in 
the semiclassical approximation. The corrections to the trace formula due to diffraction have been 
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extensively investigated in recent years. Most of the work was done within the Geometrical Theory 
of Diffraction (GTD) [^. This is an approximation in the spirit of the semiclassical approximation. 
In this approximation the Green's function is composed of free propagation from the source to the 
diffraction point multiplied by an (angle dependent) diffraction constant and followed by propaga- 
tion to the end point. Using the GTD approximation the contributions of diffracting orbits were 
calculated in a number of papers [^^^ . These contributions as well as the Berry- Tabor trace 
formula for the rectangular billiard are derived in this Appendix, for completeness. 

The GTD approximation fails in a number of cases, typically when a diffracting orbit is close to a 
classically allowed one. In these cases the contributions of the periodic orbits are given by uniform 
approximations which are more complicated, especially in the case of multiple diffraction. Such 



contributions have been examined for the penumbra diffraction for wedge diffraction ||53], for 
the diffraction from a flux line p8| and for multiple diffraction from wedges or flux lines p9|. Since 
the scattering from a point-like perturbation is accurately described using the GTD approximation, 
uniform approximations are not needed in this work. 

A convenient starting point is the Boundary Integral Method (BIM). In the case of Dirichlet 
boundary conditions this is a Fredholm equation of the second kind for the normal derivative (with 
respect to the boundary of the billiard) of the wave function. The kernel is the normal derivative of 
the Green's function of the problem with some arbitrary boundary conditions. That is. 



u(s) = -2<f ds' dnG{k; r, r')u{s') 

JdV 



(Al) 



where r(s) is a parameterization of the boundary of the billiard by arc length, dn is the derivative 
in the (outward) normal direction of the boundary, u{s) denotes the normal derivative of the eigen- 
functions and the integral is taken over the boundary of the billiard. The Green's function of the 
system is G{k;r,r'). The units h = 1, 2m = 1 are used. For the derivation of this equation and 
some applications see [^,^,^,0,^,^,0]. This integral equation has nontrivial solutions only if 



detii - Qih)) = 0, 



(A2) 



where Q{k) is the integral operator which gives the right hand side of (|A1|) when it is applied. The 
solutions ki, of this equation are the exact eigenvalues of the problem. The oscillatory part of the 



density of states is thus given by [|61|,^ 
1^ ^1 d 
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In the semiclassical limit these integrals can be approximated using the stationary phase approxi- 
mation as long as the boundary is smooth. 

The system that is investigated in this work is a rectangular billiard with a localized perturbation. 
The Green's function for the localized perturbation without the boundary is approximately given by 

Gik; r, r') ^ Go{k- r, r') + Go(fc; r, ro)D{9, 9')Go{k- ro, r'), (A4) 

where 6 and 6' are the angles of the vectors r — ro and tq — r' respectively. This approximation is 
exactly the GTD approximation and is valid far from the perturbation. The free Green's function in 
two dimensions is Go{k; r, r') = —^HQ^\k\r — r'|) ( Hq^\x) is the Hankel function of the first kind). 
The Green's function ( [A4D together with the expansion ( [A 3D lead to the periodic and diffracting 
orbits contributions to the density of states. In this Appendix it will be calculated to the third order 
in D. 

The calculation of the boundary integrals is vastly simplified by the composition law. Its semi- 
classical version is given by 

(-2)^^ / ds,--- dsn dn„Go{k; r, r„,) ■ ■ ■ d„.,Go{k; v^, n)Go{k; r^, r') ~ GI^^^; r, r') (A5) 

where the function G^J^{k] r, r') is the part of the semiclassical Green's function which results from 
the neighborhood of all classical trajectories with n bounces from the walls between r and r'. Semi- 
classically the integrals are evaluated using the stationary phase approximation and for this the 
boundary is assumed to be smooth and all the Green's function are taken in the semiclassical ap- 
proximation (for a clear discussion see [0])- Here only the simple case of one reflection from a 
straight edge will be treated to show how this method works. Consider the integral 

/oo 
dsGo{k;s,v')dnGo{k;v,s) (A6) 
-oo 

where s = (s, 0) is a parameterization of the boundary. According to the composition law 
this integral should give the Green's function that describes propagation from r' to r with one 
bounce on the boundary. We will compute this integral using the stationary phase approxima- 
tion and show that this is indeed the case. The free Green's function is given by the Han- 



kel function of the first kind and therefore (jo(^;s, r') = —jHq [kJ {x' — + y ] and also 



dnGo{k;r,s) = jkcos{a{s))H[^^ {k\/{x — s)^ + y^) where a{s) is the angle between the trajectory 



from s to r and the normal to the boundary as is shown in Fig. 
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FIG. 10. A trajectory from r' = {x',y') to r = {x,y) with one bounce 



Thus the integral to be evaluated is given by 

I{k; r, r') = H dsH^^^ {kl\s)) H^^^ {kl{s)) k cos(a(s)) (A7) 

8 J~oo 



where l{s) = y (x — s)^ + and l'{s) = y {x' — s)^ + y''^ are the lengths of the orbits from r and 
r' to the reflection point s. Since the arguments of the Hankel functions in ( [A^ ) are large in the 
semiclassical limit the Hankel functions can be replaced by their asymptotic value for large argument 



Hl^)^z) ~ J Ae*{--|--i-). (A8) 

V 7CZ 



The resulting integral is given by 



I{k; r, r') ^ ^ f rf.^^±Me^M'W+''W). (A9) 



In the semiclassical limit A; — > oo the phases in the exponent oscillate rapidly and the dominant 
contribution is from the region in which the phase is stationary, that is, near the point s which 
satisfies 

/ ^ ^ ^ • / A -1 /^\ 

sma = = = sma. (AlU) 

l'[s) l[s) 

The dominant contribution is from the neighborhood of the classical orbit which performs specular 
reflection from the hard wall. The slowly varying functions in the integral ( |A9D are replaced by their 
values at the stationary point (/ = l{s) and similarly for other values). The phase in the exponent 
is expanded to second order in the distance from s to give l{s) + l'{s) ~ Z + /' + ££li^Ii±0(^2 -yyj^ere 
= s — s. The resulting Gaussian integral can be easily evaluated leading to 

I(k- r, r') ~ , ~^ /All) 
l^kil + V) 
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This is just the semiclassical hmit of minus of the Green's function from the image of r' with respect 
to reflection around the boundary to r. For the case of a half plane this is a semiclassical version of 



the known exact result |65]. Since only the neighborhood of the classical reflection point contributes 
in the semiclassical limit it is not important that the wall of the rectangle is not infinite. Multiple 
reflections are handled by using the same method over and over. We refer the reader to where 
the case of curved boundary is also treated, for details. 

So far the boundary was assumed to be smooth, however the boundary of the rectangular billiard 
has corners. In the general case corners can lead to diffraction contributions but when the angle 
of the wedge is of the form ^ (with q integer) the only contribution is given by using the method 



of images [^]. In the case when the classical trajectory hits the corner, the resulting semiclassical 
Green's function (after integration of ([A5|) ) is the free Green's function from the image of r' under 
reflections with respect to both walls, to r. This is reasonable when one considers the method of 
images but for the angle | one can also show that it is the semiclassical limit of the double integral 
over the boundary even when the classical trajectory is reflected from the corner (the sign of the 
Green's function is positive since two reflections are involved). 

Using the composition law ( |A5|) it is now possible to compute the contributions of periodic and 
diffracting orbits to the density of states. 



1. Periodic orbits contributions 



Periodic orbits contributions are obtained when the contributions of D in the Green's function ( |A4| ) 
are ignored when substituted in ( |A3|) . Periodic orbits can be identified by two integer numbers 
(A^, M ) which denote half of the number of bounces from the horizontal and vertical sides. Therefore 
the contribution from a specific periodic orbit results from the boundary integral with n = 2N + 2M 
bounces. This contribution is thus 

^n2iv+2M-i^o(^; r2Ar+2Af , r2Ar+2M-l) ' ' ' <9n2G^o(^; r2)9„^G'o(^; 1*2, Fi) . (A12) 

The integral ( [A12| ) contains all the contributions of periodic orbits with 2N + 2M reflections from 
the walls. 

For the calculation we choose a point Si on the boundary, assume, without loss of generality, that 
it is on the lower side of the rectangle. The composition law can be used to compute the integrals 
over S2, ■ ■ ■ , S2N+2M- The contribution of the stationary phase points which belong to the orbit of 
interest is given by the Green's function that describes propagation from the image of Si under 2N 
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and 2M — 1 reflections to si as is seen in Fig. 
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FIG. 11. The stationary phase contributions for periodic orbits 



The variable Si is such that it corresponds to the point A in Fig. |Tl|. The integral over the rest 
of the variables will give contributions corresponding to the orbit of interest when S2 corresponds 
to either of the two points marked by B. These contributions are contributions of the orbit and its 
time reversal, and when either = or M = there will be only one contributing point. From the 
composition law the result of the integral over S2, ■ ■ ■ , s„ (in the semiclassical limit) is the sum of the 
Green's functions from the image points E and F to A. These contributions (and also their normal 
derivatives) are equal. For every point on the lower side there is a stationary phase contribution 
corresponding to an orbit (unless = 0). This demonstrates that the periodic orbits are arranged 
in families. The contribution to the Si integral (in ( |A12| )) from the lower side is thus 

where Ip = {2Na,j)'^ + {2MayY is the length of the periodic orbit, and (A^, M) is denoted by p. 
The factor of —4 is present since there was a factor of —2 in front of the integral (|A12|) that was not 
used in the composition law and also there are two equal contributions from the orbit and its time 



reversed one. In the semiclassical limit the dominant contribution of the partial derivative in (|A12|) 
is when it is applied to the exponent since it contributes a factor of k. The fact that dn^lp = ^M^hL 
was also used. The integrand is constant and thus the integral is just a multiplication by the side 



44 



length ttx- 

The contribution from the top side is equal to the one from the bottom side, and on the vertical 
side a similar calculation can be done, and one obtains the same contribution with the roles of 
and M interchanged. The leading contribution (in the semiclassical limit) is obtained when the 
derivative ^ is applied to the exponential in ( |A12|) . The resulting contribution of the periodic orbit 
p is thus 

di':^''\k) = 4r f ^-^^^^'"^^ ^(2M + 2iV)e-(^'^-^)) ^ ^^^^a^.m) ^^^(^^ _ ![) ^^14) 



where j3{N,M) is a counting factor taking the values 1 and |. Since the orbits with a vanishing 
index are their own time reversals, there is only one contribution for the S2, ■ ■ ■ , s„ integrals, therefore 
P(N,M) is equal to | when = or M = and is otherwise 1. 

orbits is thus given by 



The density of states measured in units of energy is d{E) = The contribution of a periodic 



(A15) 

This is exactly the contribution that can be obtained using the Berry- Tabor formula for integrable 
systems . 



2. Diffracting orbits 

The contribution of diffracting orbits, with one diffraction, to the integrals ( |A3| ) results from terms 
where for one Green's function the contribution of the second terms in ( |A4] ) is used while for the 
other n — 1 terms the first term in (|A4D is used. There are n such terms, which are all equal (they 



only differ by the labels of the integration variables). To compute the contribution of diffracting 
orbits one has to evaluate the integral 

Ji = {-2Y <j>dsids2 ■ ■ -dsndnfioik; r2,ri) ■ ■ ■ dn„_jGo{k] r„, r„_i)(9„„G'o(/c; ro,rn)D{9i,9n)Go{k; ri,ro). 

(A16) 

It will be assumed that the diffraction constant D is slowly varying so that in the region which gives 
the stationary phase contribution it can be approximated by a constant. This means that D{6,6') 
does not change considerably when the angles change by 56* ~ ./-p- where l± is the distance from 

V ^ J- 

the scatterer to the hard wall. When D is indeed slowly varying it can be replaced by its value at 
the stationary points of the boundary integrals and the composition law can be used to replace the 
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integral ( [A16| ) by the sum of all Green's functions from images of the scatterer (which cross n walls) 
to the scatterer. 

One can denote the once-diffracting orbits by the number of times they are reflected from the 
vertical and horizontal sides of the rectangle. For a pair of indices (iVi, Mi) there are four (two when 
one of the indices vanishes) such orbits as can be seen from Fig. 




The dots in the figure denote the scatterer and its images, the orbits shown are all (1, 2) diffracting 
orbits. The contributions of these orbits will result from the integral with 3 boundary points (there 
are also other contributions with A^^i + Mi = 3). Typically, there are four such orbits, but when one 
of the indices vanishes there are only two. We will add another index which runs between 1 and 4 
and counts these orbits (alternatively it is possible to allow the indices Mi, A'^i to be negative). The 
direction in which the diffracting orbit returns to the scatterer depends upon the outgoing direction 
and upon the parity of (A'^i, Mi). By counting the number of reflections on the boundary one finds 

even — even 

Oin = &{Ni,Mi)i^out) = { T . (A17) 




'(N^Min^out) -]2n- d^,,t even - odd 

odd — odd. 

In the example of Fig. |T2| the parity of (A'^i, Mi) is odd-even. The lengths of the orbits and the angles 
can be calculated easily. If the scatterer is at (a^..^^, a^^^) from the center of the rectangle then the 
lengths of odd-odd orbits are given by = ^al{Ni ± 25^)^ + a^(Mi ± 25^)^. The angles can 



46 



be computed similarly, since tan(^) = ^^l^^a^2a^5'' (^^^^ ^^at this angle is in the third quadrant if 
both the numerator and denominator are negative). The lengths and angles of other orbits can be 
computed too, just by replacing 6x by when A^^i is even and 6y by for even Mi. 

The contribution to the integral ( [A16| ) which corresponds to the {Ni,Mi) once diffracting orbits 
is thus 

jiNuM.) ^ Y.h^f'^'''D{9^u,, e(gN..,)) , ^ e-^^'-i^^'-i-) (A18) 



where /i runs over the four contributions (two contribution when either Mi or A^i vanish) and 
Ni = (A^^i, Ml) is the pair of indices that define the orbits. In what follows some of the subscripts 
Nj will be replaced by j, that is = Ij, = Oj, etc. The contribution of each one of these orbits 
to the density of states is 

J(^.MO,.(^) = ^^4: ({-ir^^''W{9i,,,Q{9i,,))^J— (A19) 



In the semiclassical limit the main contribution is from the differentiation of the exponent. Thus 
the contribution of a once diffracting orbit is given by 



JCA^i,AfO,M(fc) = (_i)A'i+M, fD(ei,,, e(0i,^))e*('='^--t'^)) . (A20) 



The contribution to d{E) is thus 



JW.M,),.(^) = (-1)^-+^^^ ^ y;^ ^ {D{9i,,, e(^i,,))e-(^'^.^-i-)) . (A21) 

The contribution of orbits with two diffractions is also needed in this work. The computation is 
quite similar, since these orbits are just a combination of two once diffracting orbits they can be 
identified by two pairs of integer indices Ni = (A^i,Mi) and N2 = [M^i^'i) and by two indices \i 
[v respectively) which denote along with the Ni (N2 respectively) orbits that are parts of the twice 
diffracting orbit. The number of boundary reflections of each segment of the orbit is denoted by 
Mi = Ni + Mi. The contribution of such an orbit to the density of states results from terms with 
two diffraction coefficients in the integral (|A3|) with n = Afi + A/2. Not all of these terms describe 
a desired orbit, since only terms in which the number of reflections between diffractions is A/i (and 
also A/2 since the orbits is closed) can describe the specified orbit. There are ^ such contributions 
if A/i = A/2 and n contributions otherwise. If A/i 7^ M there are n equal contributions that can 
be represented by an integral in which the first diffraction is, for instance, between ri and r2 and 
the second one is after A/i reflections from the boundary (the orbit is closed after additional M 
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reflections from the boundary). One can interchange A/i and N'2 without changing the resulting 
contribution. It is convenient to spht this contribution into two equal integrals in which the roles of 
Ml and N'2 are interchanged. Each of these integrals, with the diffractions between ri and Y2 and 
between r^^+i and should be multiplied by | to take into account all permutations. Therefore 

all of the contributions of orbits with two diffractions are given by 

X ^nJ^^_^,^Go{ro,r^f^+l)^nJ^^GQ{r^f^+l,r^fJ ■ ■ ■ S^jGolrs, r2)G'o(r2, ro)L'9„iG'o(ro, ri). (A22) 

The dependence of D on angles and of Gq on k was omitted in the notations for the sake of brevity. 
The sum over n is trivial. The terms of this sum (which are products of two integrals) can be 
calculated using the composition rule. The result is a sum over all possible classical paths from the 
scatterer to itself with A/i reflections for the flrst segment (and A/2 for the second). The contribution 
that corresponds to all of the twice diffracting orbits is thus 

d^oJeik)= E -^^>(-l)^^^''^^^^^'^^^(V,0(^2,.))x 

Die2,u, 0(^1,^)) ] e^(^ai.M+'2,.)-f^) (A23) 

l6lTkJlij2,u 



The dependence of the diffraction constants on the angles results from the fact that the in-going 
direction before the second diffraction is determined by the out-going direction from the flrst diffrac- 
tion and vice versa. The relation between the out-going and in-going direction is given by ([A17|) . In 
the leading order the derivative in ( |A22| ) operates only on the exponential, leading to the (energy 
dependent) density 



X 



^ (0(61,^, Q{92,,))D{e2,., e(ei,^))e^('=('^-+'^-)-t'^)) . (A24) 



Combining ( |A24 ) with ( |A15 ) and ( A21 ) we flnd that the total density of states taking into account 
terms up to the second order in D is given by 



d^^m = E ^^^^^^e^^'^^-f ) + E i-if^^''^-^^D{e,,„ e(0i,,))e^(^'^-i-) 

N nJSirklp Ni,/. AirkVSTck 



+ c.c. (A25) 



+ E E (-1)^1+^^!+^^^^^ ^''^ +;^'- fD{9i^„ e(02,.))/^(^2,., 0(^^i,,))e^('^('^-+'^-)-i-)) 
Ni,/.N2,!/ G4:n'^k^Jh^f,l2,u 
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Note that when one of the indices of Ni (N2 respectively) vanishes the sum over /i (z/ respectively) 
includes two contributions instead of four. Thus for an angle independent diffraction constant such 
orbits will have a factor of half similar to the one that appears in the periodic orbit contributions. 
Equation ( [A25| ) will be used to compute the various contributions to the spectral from factor. 



Similarly one can calculate the contributions of orbits with three diffraction. The calculation 
follows the one of two diffractions closely, and the only difference is that one now has a factor of | 
instead of | in ( [A22D . Since there are three diffractions the sum is over Ni,yUi, N2,/i2 and 1^3, fis 
which describe the three segments of the orbits. In this work the contribution to third order in D is 
required only for angle independent scattering. In this case the contribution of all orbits with three 
(angle independent) diffractions is found to be 



(^E) = ^ ^ y^^ ^_i-^Ni+Mi+N2+M2+N3+M3 ^ h,^ll + h,^J.2 + h,tA.3 ^ 
Ni.^ti N2,At2 N3,At3 127r/c(87rA;)2 Jll,|^J2,^l2h,tM3 



X ^iHh.n+h.^.2+h,n) + c.C. (A26) 

In Sec. it is assumed that the scattering is angle independent, therefore the function D{6,6') 
is replaced by an angle independent constant D. Moreover we are interested in the form factor 
where in the calculation energy averages are taken and sums over periodic orbits are approximated 
by integrals. Within such approximations the weight of the orbits with I3i^n,m) (*-"^® vanishing 
index) is negligible, therefore in other parts of the paper these factors were replaced by unity. Within 
these assumptions (|A25|) and (|A26|) reduce to (H). 



APPENDIX B: THE DIFFRACTION CONSTANT OF THE POINT SCATTERER WITH THE FREE 

GREEN'S FUNCTION Go 

The diffraction constant is the on shell transition matrix, T, which for a point interaction is given 
by 

z — iK /■ ,5 ^ , „x , z + iA. 

1 - e*«i 

where ^y/z > and Go is the Green's function in free space. To compute the integrals we use the 



D-Hz) 



■ / d^x Go(z; X, 0)Go(iA; x, 0) + ^ ^ / (fx Go(z; x, 0)Go(-iA; x, 0) (Bl) 
J 1 — e~*" J 



following integral representation for the free Green's function in two dimensions 

G„(.;x,x„)=p-^/d=*^-^. (B2) 



Let us compute the integral 

£x Go{z; X, 0)Go(^A; X, 0) = j—:^ J Ski d'^h £'x _ .^^ 



iki-x g— jk2-x 



TT, TIT, TT- (B3) 



(27r)2 7 (F -z)(A;2 -zA)' 
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Changing variables to E = k"^ one obtains 



1 

I = — I dE — = — ; — / dE 

^ in Jo {E-z){E-iK) ATx{z-iK)Jo 



E-z E-iA 



Using the known relation 



P- 



1 



=F iiT6{x — xq) 



X — xq ± ie X — xo 

and noting that we are interested in the case where Qz tends to zero from above we obtain 



1 



lim lim 



2-/3 I j-M \ j-M 

dE + / dE / dE 

E — z Jz+p E — z Jo 



Att{z - iA) 

The integrals are elementary, and taking the limits results in 

1 





tt' 











Att{z - iA) 

The second integral in ( pi]) is computed in the same way 

1 



J d^x Go{z; X, 0)Goi-iA; x, 0) 



4:7r{z + iA) 



In 



Substituting ( P7| ) and ( P8|) into (|B 1| ) leads to the diffraction constant: 

D{z)- 



2 ^'Ha; 4 1-cosa ^ 2 



E-iA 



+ in 



.Stt 



(B4) 



(B5) 



(B6) 



(B7) 



(B8) 



(B9) 



APPENDIX C: ASYMPTOTICS OF SCATTERING FROM LOCALIZED POTENTIALS 

In this Appendix Dg{q) and hCk, q) are computed up to order k'^a'^. The computation is presented 
explicitly only to order ka. The calculation of the second order is straight forward and fairly similar 
to the calculation of the first order, but its details are too long to be included in the paper. 



1. Calculation of ^^^(q) to the order (ga)^ 

In the calculation of Dg{q) many of the contributions cancel, but equation (^) does not show 
this explicitly. A more convenient way to compute -Dg(q) is obtained when (^) is iterated, and all 
the resolutions of the identity in momentum space which appear there are replaced by resolutions 



of identity in coordinate space. Using (|83D leads to 



Vo 



(a) 
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e^'^q-ya (CI) 





i-l n \ f J2 (a)/(yi)/(y2''^) iaq-(yi-y'"')r^ /; (a) \ 1 /(Y: 

X U{yi)Go{ka;yi,y2) - J d yy — e "^^^^ ^Go(A;a; ^ ya) I — ^ 

X (|/(yi)Go(fca;y,,y,) -/ /(yi)/(y^"\ ...,y.^yr')g,(fc,; yW, 

X (^/(y.)Go(^a;y„y3) -|^^,(-) ^^y-)^Jy3°X --c,,y.-yr')g^(,,. y(-), j /i| 
where 

K, = (q|U|q)= Jd'yfiy). (C2) 
The terms of higher order in g have the same structure, but with more of the factors such as 

B, = f{y,)G,{ka- y„ y,+0 - / rf2^(a) /(y.)/(yj+\) ^.aq.(y,-y(;),)^^(^^. y y^.^^). (C3) 

Equation (|Cl| ) is a convenient starting point from which the asymptotics of Dg{(\) in the hmit qa ^ Q 
can be calculated . 

To calculate the leading order of -Dg(q) (with respect to qa) all the exponentials are replaced by 
unity and the Green's function by the leading order as given by (PO]). The leading order of any of 
the brackets is given by 

B, - /(y,)0(y„ y,+i) - / d^fl J^^'^^^^'^^'\ {y% y,+0- (C4) 

Note that the contribution from the constant parts of the Green's functions, di, cancel each other. 
This is the advantage of the Noyce method when describing scattering at low energies in two di- 
mensions. Most of the Green's functions in equation ( |G1D appear in brackets of the form Bj, and 
thus, in the leading order, their constant parts cancel. However, in every term of some power in g 
there is one Green's function that does not appear in some brackets. When this Green's function is 
multiplied by one of the brackets one finds that the contribution of its constant part, cii, vanishes. 
This is true since for any functions 5'i(y),fi'2(y)5 G(yi,y2) one finds 

/ d^yt^d^^d^, g,iy[^^)d, ||/(yOG(yi, y^) - / d^^^ l^y^^^I^G{y^^\ y,)^ g,{y,) = 

J d^y^^d'y, (7i(y!"Vi (/ Ai/(yi)G(yi, y^) - / d^y^^ f{y^^^)G{yt\y2)) g2{y2) = 0. (C5) 

Therefore the constant part of the Green's function does not contribute to the leading order of -Dg(q) 
except for the term of order g. Using ( p^) and (|C5|) the leading order (in qa) contribution to Dg{q) 
is given by 
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/(yi 

Vn 



-g^Jd'yi'^'^d'yA /(yJ'^^)0(yJ"\ Yi) /(yi)0(yi, y2) - /ci^i/^' 



2 (a)/(yi)/(y2"^) (a) /(y2j 



-0(y2 ,y2) 



-g'Jd'yi^^d'y.d'y^d^s /(yi"^)0(y!'^\ yi) [^/(yi)0(yi, y2) - /^'y^' 
X /(y2)0(y2,y3)- Jd'y^^ 



(a)/(y2)/(y^^) (a) /(y 
— 77 (pyy-i ,y3 



(C6) 



Note that the integrals in ( |C6| ) are angle independent. To simplify the notations the following 
definitions are useful 



Lo = J Ai /(yi) = Vo 

= j d^yid^y2 /(yi)0(yi,y2)/(y2) 
^•2 = 1 f^^i/io^Vo?^2/3 /(yi)0(yi,y2)/(y2)0(y2,y3)/(y3) 



(C7) 



L2L1 L\ 



(a 



Using these definitions (CC) can be written as 

Vo V I'D I'D / \ ''0 I'D ''0 , 

Note that the contribution of order is a sum where each term is composed of a product of Lj's 
where the sum of j's is A^. It is multiplied by — to the power of number of i^j's. The sum is over 
all possible partitions of into positive integer numbers so that X]j = Such partitions 
are generated by the expansion 

1 k 



Vig) 



Vo 



Vn 



J2m=0 g^Lm Vo + Y^m=l g'^Lm 

that can be written in terms of partitions as 

00 li 

v{g) = i+Y.9''Y.Y{ 

N=l i j=l 



1 + E 

k=l L m=l 



00 T 

Vo 



Vo 



(C9) 



(CIO) 



where i is the index that specifies the z-th partition consisting of k terms. The various indices of the 
z-th partition satisfy Z^jLi f^j = N with the requirement that all rrij are positive. Consequently the 
leading order of Dg{q) is 



Dg{q) = V{g)-gd^Vo. 
Expanding V{g) to the order one recovers 



(Cll) 
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There are no contributions to Dg{q) of order qa and the first correction to Dg{q) scales as g^a^ 
(times some power of Inga). To see this note that from ( pSD it is clear that -Dg(q) is built from 
products of matrix elements of the form (q|UGUGU ■ ■ ■ UGU|q) with various numbers of Green's 
functions. A contribution of the order qa must come from one of these elements which is expanded 
to that order. For example, the contribution of (q|UGU|q) to this order is 

-J (fyi(fy2 iaq- y if {yi)Goika;yi,y2)f {y2) + J ^V^^V ^aq ■ y2/(yi)G'o(fca; yi, y2)/(y2) =0, 

(C12) 

since Go{ka; yi, y2) = Go{ka; y2, yi). This holds if more Green's functions are used, therefore Dg{q) 
cannot have contributions of the order qa. 

The contribution of higher order in qa to -Dg(q) can be computed by expanding the exponents and 
the Green's functions in (|CT| ) and collecting all the terms up to the required order of qa. Here -Dg(q) 
is calculated to order q'^a^. The calculation is straightforward but tedious and only the final result 
is presented 



Dg ^ Vig) - gd.Vo + (d,-^] VoVig) 



oo r 

E - Vig) 



fc=0 ^0 



1=1 Vo' 



oo 2i0 (++) I (+-) , (-+) , -2ie {—) 
-q'a' (Vig) - gd.Vo) (Vig) - gd.Vo) E ^' ^' ^' 



1=1 



4K) 



9 



gdi 
2 



Oq + -e (Jo + -e (Tq 



+ V{g)gdiVo E —r -9'' 



1=1 

, 2 



Ol +2^ + 9^ 



l„-2i 



2Vo 



"9^ 



\k=0 



2U, 



(C13) 



where 9 is the direction of the vector q in the plane. The constants Ck, cik, ^k, CTfc and /i^"''-* are 
independent of 6 and q and are defined in what follows. First we define 



^k = J d^yi ■ ■ ■ d^Vk+i (riiyj - yi+i)^j /(yi)0(yi, y2)/(y2 



) X ■ • • X 0(yfc, yfe+i)/(yfc+i) (C14) 



and 



ao= I Aie"^'^|yi|/(yi) 
«i ^ y dVA2 e-*^^|yi|/(yi)0(yi,y2)/(y2) 
^2= c?V'^^2/2rf^2/3 e"''^'|yi|/(yi)0(yi,y2)/(y2)0(y2,y3)/(y3 



(C15) 
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where 9i denotes the direction of yi. The other constants are 



5, 



/ Ai |yir/(yi) 
5i = j d?yid^y2 |yi|V(yi)0(yi,y2)/(y2) 

h = j d^yid^y2d^yz Iyir/(yi)</'(yi,y2)/(y2)(/'(y2,y3)/(y3) (C16) 



a, = j dVy2 e-'^^^|yi|V(yi)0(yi,y2)/(y2) 

a^^J d'y.d'y^d'y, e"^^''^ |yi| V(yi)0(yi, y2)/(y2)0(y2, y3)/(y3) (C17) 



and 



(++) 



I d'y^d'y, e-^^^|yi|/(yO0(yi,y2)/(y2)|y2|e-^^^ 
^ I AiA2e-^^^|yi|/(yi)0(yi,y2)/(y2)|y2|e+^^^ 

^ I A1A2A3 e-^''Myi|/(yi)0(yi,y2)/(y2)0(y2,y3)/(y3)|y3|e-^^^ (C18) 
4--) = J d'y^d'y^d'y^ e^'^|yi|/(yi)0(yi,y2)/(y2)0(y2,y3)/(y3)|y3|e^'^ 

where the lower index specifies the number of 0's in the integrals and the upper index is the opposite 
of the sign of the angles in the exponents. Note that 

/^r^=^r^* , f^-'=^''* (C19) 

In addition to these definitions, for future use, it is also convenient to define the sums 

00 



00 

fc=o 

00 



K9) = E ^fc^' 

fe=0 

00 

o{g) = Y.o^g'^ (C20) 

fc=0 

00 



fc=0 
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2. Expansion of h(k, q) for ka and ga ^ 1 



The most convenient way to obtain the expansion of h(k, q) for small ka and qa is by iterating 
and inserting resolutions of the identity (in coordinate representation) where appropriate. Note that 
by definition (^2|), /i(q, q) = 1. Changing the integration variables to yj = ^ results in: 

/i(k,q) = J d^y,e- -^-i^~^yy-^-^ + g J d'yod'y, |^e-*"'^-y«/(yo)Go(A;a; yo, yi) 
J Vq / Vo 



+ ■■• (C21) 
The general term in this series is composed of a first factor 

^0(y0,yi) = e-^'^''-^V(yo)Go(fca;yo,yi) - / rf^^W /(yo)/(yS"^) ^.a(q-k).yog^.aq.y<")^^(^^. 



/(y2 



X |/(yi)Go(A;a;yi,y2)- / rf^^f^:iI£±i^^A:i^e^'^''-^^-^^ )G,{ka-y'^\y,) \ ^e^'^'^-^^ 



(C22) 



followed by several factors of the type Bj defined by (|C^ ), and a last factor 

= e^'i-y^^^^^ (C23) 

(with all variables integrated). To compute the ka and qa dependence of /i(k, q) every term in 
equation (|C21| ) is expanded for small ka (qa) and the whole series in g is summed. In what follows 
/i(k, q) will be computed to order /c^a^ (and also g^a^, kqa^). 

To compute /;.(k, q) for small ka the first factor has to be calculated at least to the first order ka. 
In this order 

^(yi) = j c?^yo^o(yo,yi) 

j Ao^a(q-k).yo (|/(yo)0(yo, yi) - / ^^^^^^^^^^^ • (C24) 

Note that in the leading (zeroth) order in ka this contribution vanishes. This is related to the 
fact that by definition /i(q, q) = 1, and in the leading order the direction of the momentum is not 
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important. To compute h{k,q) up to order ka the factors Bj and An can be taken in the leading 
order (zeroth order in ka). Using and ( |C24| ) one finds 



h(k, q) = 1 + m jy" d'^yo (q - k) ■ yo 



/(yo) 
v. 



a)/(yo)/(yi 



-0(yi ,yi 



+ ^ y" d^vod^yi (q - k) ■ yo j^/(yo)0(yo, yi) - J d^yt^ 
+ 9' J d^od'yA (q - k) ■ yo (|/(yo)0(yo, yi) - / ^^i/i"^ ^^^°^^J^^"^V (yS"\ yi 
X |^/(yi)0(yi,y2)- J d'yi^^ 



/(yi 



a)/(yi)/(y2 



-0(y2"\y2 



/(y2 



+ 



(C25) 



Since k ■ yo = A;|yo| cos(6'' — Oq) (and q ■ yo = qIyoI cos{9 — Oq)) the dependence on the incoming and 
outgoing momentum can be separated from the integrals in (|C25|) . Use of ( P15|) results in 



h{k, q) 



l + 'T ( 



ie'\ "0 



( 



Vo 
fje 



e — e 
L ^2 



^2 _ -^1 



Vo Vo 
K) Vo 




Vo V 



e - e- 



Vo \Vo VI 



Vo Vo 
/\ a 
Vo 



ie'\ ^ 
Vo 
ie'\ ^ 
Vi 



+ ■ 



(C26) 



First note that ao appears with the series 1 — — g'^ — ^ 



■ ■ ■ = V{g). This results from 
the fact that to get the ao contribution from the second term in the first brackets in ( |C25|) it must 
be taken in all orders of g^ . Then the integrals from y^'^'' to yN behave exactly as the leading order 
(in qa) contribution to Dg(q). Similar considerations can be applied to each of the aj. This series 
(in g) is summed to give 



ka 



l + ^'^V{g)^^{e^'-e^'' 



k=o ^0 



00 * 



Vo 



k=o ^0 
Vo J 



9 



(C27) 



with a{g) defined by ([C20D . 

The computation of the contribution to h(k, q) of order fc^a^ (with additional powers of In ka) 
follows the same lines. However, it is too tedious to be included here and only the result is given, 



Mk,q).l + .^n.){(e^^-e-)^ + ( 



e^'-e- 



Vo 



+ k'a' {y{g) - gd^Vo) 



CO I 

1=1 4>/o 
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2Vr 



{gVod, - V{g)) 







1 



Oq + -e (To + -e cr. 



1 







-2ie -2i6 



i9'\ ( ^*{9) 
2K) 



2K, 



2K, 



(C28) 



3. The angular dependence of D{6' , 



Finally the diffraction constant of (|78| ) is computed to the order k'^a^ with the help of (|86D . It is 
useful to expand (|C13| ) as: 



oo A 



1 + A:^a^ ^ iy{g) - gV,d,) ^ ^ (e^^VS^^^ + /^["'-^ + /^S"^^ 



.=1 4K, 



, „-2ie,,{ — )\ 9di f J- ^,,2ie„ , ^,.~2i0*\ 



(5; + ie^^V, + ie-2^v* 
2U^ 



+ V{g)gd,Vo E .T. 

\fc=0 "^*^o 



«=1 '"'0 



where Q(c/) is defined by O- Substituting (|C28D and (|C2|) in (H) leads to (|9 



(C29) 



APPENDIX D: CALCULATION OF ANGULAR DEPENDENT CONTRIBUTIONS OF THE FORM 

FACTOR 

In this Appendix the details of calculations leading to some of the results of Sec. [V^ are presented. 

1. Diagonal contributions to the form factor for a scatterer at the center 



For odd-odd orbits the incoming direction is given by 7r + where 9 is the outgoing direction ( |A17| ). 
The contributions from the odd-odd orbits in (|l^) is replaced by 



^(l.o-o)( 



^2 „i 



327r Jo 



■ de \D{9, 9 + n) + D{9 + n, 9) + D{tx - 9,2n - 9) + D{2n - 9, 



(Dl) 



r 



2 r27T 



327r 



/ "d9 D*(9, 9+7t) [D(9, 9+n) + D(9+7i, 9) + D(n-9, 2ti-9) + D(2n-9, n-i 
Jo 



The contribution from the odd-even orbits is given by 
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K(i'°-^)(r) = — f'de \D{e,TT-e) + D{n-e,e) + D{7r + e,27r-e) + D{27v-e,7v + e)f (D2) 

32tt Jo 

[^^d9D*{9, 71-9) [D{9,n-9) + D{tx-9, 9) + D{'k+9, 2'k-9) + D{27i-9, 7i+9)] . 
Jo 



T-2 r27T 



32n 

The contribution from the even-odd orbits is 



2 „IL 

j^{i,e~o) (r) = rd9 \Di9, 2n - 9) + D(2n - 9,9) + D(n + 9, n - 9) + D(n - 9, n + 9)f (D3) 
32n Jo 



327r Jo 

T-2 r27T 



/ d9D*{9, 2-K-9) [D{9, 2^-9) + D{27i-9, 9) + ^(vr+e, ■k-9) + ^(vr-e, 7r+^ 

JO 



327r 

The diagonal contribution of periodic and diffracting orbits to the form factor (up to order r^) is 
obtained by summation of (KI^) and ( pi[ - p3| ), resulting in 



7^(1) (r) = 1 + -^$5 r d9 D{9, 9) 

271 Jo 

+ — (^^ d9 {D*{9, 9) [D{9, 9) + D{7t - 9,7i - 9) + D{7t + 9,7r + 9) + D{27t -9,27t- 9)] 
32n Jo 

+ D*{9, 9 + 7i) [D{9, 9 + 7r) + D{9 + vr, 9) + D{7i -9,2^-9) + D{27i -9,7i- 9)] 
+ D*{9, 7T-9) [D{9, 71 -9) + D{7i -9,9) + D{7i + 9,27^-9) + D{27t -9,7i + 9)] 
+ D*{9, 271 - 9) [D{9, 27T-9) + D{27t -9,9) + D{7t + 9,7t-9) + D{7t -9,7t + 9)]} . (D4) 

Substitution of (H) yields ( |105|) . 

2. Non diagonal contributions to the form factor for a scatterer at the center 

For twice diffracting orbits the diffraction constants appear in the amplitude in the form 

^ D{9,,^, Q{9,^,))D{9^^,, 0(9,,,)) (D5) 

instead of that was used for angle independent scattering. The functions of ( [A17| ) depend on 
the parity of Ni and N2 but since on the saddle manifold Ni + N2 = Np, where the indices of Np, 
corresponding to a periodic orbit, are all even. Therefore both Ni and N2 have the same parity and 
for orbits that will contribute both the functions G are identical. 

The computation of the contribution of the saddle manifold can be done exactly as in Eq. (pOD -(p5|). 
Since in the semiclassical limit the angles of orbits which contribute are given by 9p + 69 where 9p is 
the direction of the periodic orbit and 69 1, the expression (P5|) can be taken out of the saddle 
integral because it changes slowly on the saddle manifold. The various non diagonal contributions 
from any type (even-even, even-odd, etc.) of twice-diffracting orbits and periodic orbits are given 
by: 
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)(r) = ^de^Yl D{e„ Q{e.))D{e^, Q{e,)), (D6) 

A TT Jo 16 ^ 

where x — x denotes the parity of Ni (and also of N2), that is, of the segments of the twice diffracting 
orbits. 

For even-even orbits, from (|A17|) one finds B(^) = 6. For each segment that starts in the direction 
9 where < < | (or fi = 1) there are also segments of identical length in the directions ix — 9 

{fi = 2), 71 + 9 = 3) and 2tt — 9 (/x = 4). Therefore there are 16 contributions to the sum in ( 

2 /.z 

K'(^^^-^)(r) = -3?— / ' {D{9,9)D{9,9) + D{9,n - 9)D{7i - 9,9) + D{9,7r + 9)D{7r + 9,9) 
32n Jo 

+ D{9, 27r - 9)D{2ti -9,9) + D{n - 9, 9)D{9, n - 9) + D{tt -9,n- 9)D{n -9,n 
+ D{n ~9,n + 9)D{n + 9,n-9) + D{tt - 9,2n - 9)D{2tt -9,tt-9) 
+ D{n + 9, 9)D{9, n + 9) + D{Tr + 9,tt - 9)D{tt -9,n + 9) 
+ D{n + 9,TT + 9)D{tt + 9,n + 9) + D{tt + 9,2n- 9)D{2tt -9,n + 9) 
+ D{2n - 9, 9)D{9, 2tt - 9) + D{2n -9,n- 9)D{n - 9,2tt - 9) 
+ D{2tt -9,tt + 9)D{n + 9,2n-9) + D{2n -9,2tx- 9)D{2tt - 9,2tt - 9)} 



r d9 {D(9, 9)D(9, 9) + D{9, vr - 9)D{n -9,9) + D{9, tt + 9)D{7t + 9, 
Jo 



327r 

+ D{9,2TT-9)D{2n-9,9)}. (D7) 

The calculation for diffracting segments with other parities is similar. For the odd-odd orbits O(^) = 
TT + 9. After a similar calculation the non diagonal contribution from these orbits to the form factor 
is found to be: 

Js:'(2-°-°)(r) = -—3? d9 {D{9, TT + 9)D{9, n + 9) + D{9, 2^ - 9)D{'k -9,71 + 9) 
32tx Jo 

+ D{9, 9)D{ti + 9,^ + 9) + D{9, tt - 9)D{2tt -9,7r + 9)}. (D8) 
For orbits of the odd-even type, Q{9) =71 — 9 and their contribution is given by 

K (2'°-^) (r) = 3? / d9 {D{9, n - 9)D{9, n - 9) + D{9, 9)D{tt -9,tt-9) 

327r Jo 

+ D{9, 2tx - 9)D{tx + 9,tx-9) + D{9, tt + 9)D{2t\ -9,tx-9)}. (D9) 
And, finally, for the even-odd orbits, O(^) = 21^ — 9 and 

^'(2,e-o)/ ^) _ -—^ d9 {D{9, 2tt - 9)D{9, 2-^-9) + D{9, 9 + n)D{n - 9,2n - 9) 

3277 JO 

+ D{9, TT - 9)D{tt + 9,2^-9) + D{9, 9)D{2tv -9,2tv-9)}. (DIO) 

The total non diagonal contribution to the form factor is obtained when (|D7| )- (P10|) are summed 
leading to 
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2n 



327r Jo 

+ D{e,n -e){D{e,TT 



3? / de [D{e, e) {D{e, e) + D{n-e,7v-e) + Din+e, n+e) + D{27v-e, 27v-e)} 
+ D{7T -e,e) + D{7v + e,27v-e) + 0(271 -e,-K + e)} 
+ D{e, -n + e) {D{e, n + e) + D{n + e,e) + D{n -e,2TT-e) + d{27t -e,n-e)} 

+ D{9, 2tt - 9) {D{e, 2TT-e) + D{2n -9,9) + D{n + 9,tt-9)+ D{n -9,n + 9)}]. (Dll) 



Substitution of (|93D leads to (|108|) 



3. Diagonal contributions to the form factor for a scatterer at a typical location 



Odd-odd orbits with different indices fi differ in length and for every /i there is a separate contri- 
bution to the diagonal part. The sum of these four contributions is given by 

^2 



r 

32n Jo 

T 



'd9 {^\D{9,9+Ti)\^ + \D{'k-9,2ti- 
d9 \D(9,9 + 7i)\\ 



+ \D{9+7c,9)\^ + \D{27i-9,7r-9)\^) 

(D12) 



327r Jo 

An orbit of the odd-even type which starts at the direction 9 will return to the scatterer at the 
direction tt — 9. Its length is identical to the length of its time reversal which starts at the direction 
271 — 9 and returns with n + 9. Thus orbits of this type with ^ = 1 and /i = 4 (with the same 
index Nj) have the same length. This is also true for orbits with n = 2 and /i = 3. These length 
degeneracies lead to a diagonal contribution to the form factor which is given by 

^2 

D{2n - 9,TT + 9)\'^ + \D{tt- 9,9) + D{n + 9, 2n 



T 

327T Jo 



d9 [\D 



TT 



d9 D*(9,n 



[D{9,n 



D{27i-9,n + 9)] 



(D13) 



327r Jo 

The contribution of the even-odd orbits is very similar to the one of the odd-even ones. The difference 
is that orbits which start at direction 9 return at the direction 2tt — 9 and their time reversal starts 
at TT — Therefore the lengths of orbits with /i = 1 and /i = 2 (and the same index Nj) are identical 
and it is also true for orbits with = 3 and = 4. Their contribution to the diagonal form factor 
is given by 

^2 

\^ + \D{2'R -9,9) + D{-n + 9,1^ 



_^(l.e-o)/ 



r 

m Jo 



d9 (\D{9.2tt 



D{ 



TT 



9, 



TT 



f d9 D*{9, 27T - 9) \D{9, 2ti - 9) + D{n - 9, 
Jo 



TT 



_ , , (D14) 

327r JO 

The diagonal contribution to the form factor is obtained when ( |104] ), ( pi2|) , ( |D13|) and ( pi4|) are 
summed, resulting in 
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f-27r 



+ — I dO {D*{e, 6) [D{e, 6) + D{-K -6,71-6) + D{-n + 6,ti + 6) + D{2ti -6,27c- 6)] 
P + D*{6, TT-6) [D{6, n-6) + D{2n -6,n + 6)] 
D*{6, 2n - 6) [D{6, 2ti-6) + D{n -6,n + 6)]}. (D15) 



1 + 


— <y 

2-n 






327r 


i 


m 




D*( 


6. 2n 



Substitution of (H) leads to (|TTOD 



4. Non diagonal contributions to the form factor for a scatterer at a typical location 

For odd-odd orbits if the first segment has /i = 1 then the second segment must have /i = 3 to 
satisfy (|56|). Similarly if the first segment has ^ = 2 then only segments with = 4 can contribute 
(and vice versa). Therefore, instead of the 16 contributions for (P8| ) when the scatterer is at the 
center of the rectangle there are only 4 contributions for a scatterer at some typical position. The 
resulting contribution to the form factor is given by 

2 2L 

ir'(2'°-°)(r) = -— d6 \D{6,6)D{6 + n,6 + 7i) + D{n - 6,71 - 6)D{2n - 6,2n - 6) 
327r Jo 

+ D{6 + 7r,6 + tt)D{6, 6) + D{2^ -6,27i- 6)D{it -6,7i- 6)] 

= -7^3? d6 D(6, 6)D(6 + tt, ^ + vr). (D16) 
327r JO 

When both segments are of the odd-even type there are length degeneracies between a segment and 
its time reversal. If the first segment has /i = 1 then the second segment can either have = 3 or 
fi = 2 leading to the combinations of the form D{6, 6)D{ti -6,71-6) + D{6, 2-n - 6)D{ti + 6,71-6). 
Similar contributions are obtained when the first segment has fi = 2,3,4. Thus instead of the 16 



contributions leading to (|D9| ) when the scatterer is at the center there are only 8 when the scatterer 
is at a typical position. These contributions are given by 

K'(^'°-^\t) = 3? / ' d6 [D{6, 6)D{ti -6,n-6)+ D{6, 2tt - 6)D{tt + 6,n-6) 

32tc jo 

+ Di-K -6,71- 6)D{6, 6) + D{tt - 0, tt + 6)D{2tt - 6, 6) 

+ D{n + 6,TT + 6)D{2n -6,2tt-6) + D{n + 6,n- 6)D{6, 2tt - 6) 

+ D{2n -6,2n- 6)D{n + 6,tt + 6) + D{2tt - 6, 6)D{tt -6,tx + 6)] 

= -—3? d6 (D(6, 6)D(n -6,n-6) + D(6, 2n - 6)D(7i + 6,n-6)). (D17) 
327r Jo 

The contribution of the even-odd twice diffracting orbits can be computed in a similar manner. Here 
if the first segment has yU = 1 the second is either /i = 3 or yU = 4, with the same length, resulting in 
the non diagonal contribution to the form factor from these orbits. 
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K (2.^-")(r) = K / de {D{9, 9)D{27i - e,27r - 6) + D{9, n - 9)D{n + 9,271- 9)) . (D18) 

3277 JO 

The non diagonal contribution is obtained by summing (p7|) and ( pl6| - |DT8| ) resulting in 



K'(2)(r) = -— 3? r d9 {D(9, 271 - 9) [DU + 9,7i-9) + D(27v - 9, 9)] 
327r Jo 

+ D{9, 9) [D{9, 9) + D{7T -9,7i-9) + D{7i + 9,7i + 9)+ D{27i -9,2'k- 9)] 

+ D{9, 71-9) [D{7i + 9,271-9) + D{7v - 9, 9)] + D{9, vr + 9)D{7t + 9,9)}. (D19) 

Substitution of (|3|) yields (|IT^). 

APPENDIX E: COMPARISON TO THE PAPER OF BOGOMOLNY AND GIRAUD ^ 

In this Appendix results of the paper by Bogomolny and Giraud are compared with the results 



of our present work where appropriate. The notations in this Appendix are the ones used in |^ . 

Consider first the case of angle independent scattering. For periodic boundary conditions the result 
of is given by (61) and (62) there. Expanding to order and using 



Aoo2 = 1 (El) 

leads to 



1 1 /I 1 , 



K^,(t) = 1 - -\D\'t + —iDl^r' + — IDI* - — l^r t\ (E2) 
pbcK J 4' ' 32' ' V32' ' 384' ' / ^ ' 

The form factor for periodic boundary conditions is related to the form factor for (an angle inde- 
pendent) scatterer at the center by the relation ([T5|) of our paper. Substituting (^) for Kper in ( ffSD 
leads to the result obtained in this work for a scatterer at the center (|52|). 

When the scatterer is at a typical position the results are given by Eqs. (83), (84) and (85) of |p8|| . 
Expanding to order one finds 

K{t) = 1 - \\D\'{r)r + ^\D\'{r'y - ^\D\'{r'y + ^A,o2\D\'t\ (E3) 

Using 

(r) = l 
^ ' 4 

ir') = I (E4) 

and 



62 



81 

^02 = - (E5) 
lb 

leads to an expression for the form factor that is identical to the one computed in our work and is 
given by (p7|). 

The form factor of an angle dependent scatterer at a typical position is given by Eqs. (109) and 
(110) in It is sufficient to expand it to order for comparison with the form factor computed 
in our work. Eq. (109) in leads to 

K{t) = ((1 + \-SS{cp, cp)r + 1|5(^, ^)| V - j^^S'i^, ^y)) + Oir% (E6) 

where 

S{ip, ip) = D{(p, (p) + D{tc - (p,TC - (f) + D{-n + (p, -TT + (p) + D{-(p, -cp) 

- [D{ip, + D{n - V5, -TT + ^)] e^'^' - [D{-ip, ip) + D{-n + ip,n- ip)] e'^'^' 

- [D{^, n-^) + D(-v9, -TT + ^)] e'* - [D{-n + ^, + D{n - ^, ^)] e-^'^' 
+ D{ip, -vr + yp)e2^('^^+'^^) + D{-ip, vr - ip)e^'^'^'-'^^'^ + D{-7i + (p, !p)e-^'^^'+^^'> 

+ Din - ip, -y,)e2^(^2-<^i) (E7) 

is obtained with the help of (88), (102) and (105) there. The brackets ((■ ■ •)) denote averaging over 
the uniformly distributed random variables 0i,02 (see the paragraph following (79) of |^), and 
over the angle ip in the range (0, |). The averaging over 0j will leave only the terms without the 
exponentials. We show that the terms in (^) are identical to the contributions to the form factor 
computed in our work. For instance, the term linear in r is proportional to 

{{Siip, ip))) = - r dip [D{ip, ip) + D{TT-ip,TT-ip)+ D{-Tl + ip,-TX + ip)+ D{-ip, -ip)] (E8) 
TT Jo 

which is identical to the contribution linear in r in ( pi5|) of our paper (see also (|102|) ). The term 
resulting from jS'p is explicitly given by 
1 1 f~ 

— {{\Siip, ip)\^)) = -— "dip (\D{ip, ip) + D{n~ip,Tr-ip) + D{-n + cp,-n + ip) + D{-ip, -ip)\^ 
64 Sin Jo ^ 

+ \D{ip, -ip) + D{n -ip,-n + ip)\'^ + \D{-ip, ip) + D{-tx + ^, tt - y?)p 

+ \D{ip, TT - V^) + D{-^, -TT + <^)P + \D{-'K + V5, -^) + D{'K - (^)p 

+ |Z}(<^,-7r + (^)p + |D(-yp,7r-yp)|2 + |Z}(-7r + <^,<^)P + |Z}(7r-<^,-V9)|2) (E9) 



and is identical to the diagonal contribution from once diffracting orbits, that is, the prefactor of 



in (D15) as can be easily be seen examining the preceding equations. Finally 
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' d(p {[D{(p, if) + D{Ti-ip,Ti-Lp) + £)(-7r + V5, -vr + v?) + D{-ip, -9?)]' 
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+ 2 -^) + Z}(7r - V?, -TT + V?)] [/^(-^, + D{-n + ip,n- ip)] 

+ 2 [D{if, Tc -if) + D{-ip, -TT + ip)] [D{-7T + ^, + D(7r - V?, 
+ 2D{ip, -TT + (^)D(-7r + v?) + 2D{-yD, n - p))D{'n - v?, -^) 



(ElO) 

After some manipulations ( |li)lCI| ) can be shown be identical to ( pi9| ) of our paper. 

If the scatterer is not at a typical position the calculation can be done within the formalism of 
P^ , but the average over 0i and 02 has to be modified In particular for the scatterer at the 
center the form factor is obtained from ([tL)3|) and (|bLl6|) when the average over 01, 02 (defined by (71) 



and (72) of [^) of a function / is defined by 

1 



(/(01,02)) 



E E /('^i'</'2). 



(Ell) 



For an angle independent scatterer, with the help of (83) and (85) of it leads to 



where n 



1, 2, ■ ■ ■ and to 



A 



002 



16. 



(E12) 



(E13) 



Substitution in (|E3|) results in the form factor for a scatterer at the center (|52|) . 

The form factor for an angle dependent scatterer at the center can also be computed applying the 
average ( [Ell| ) to (|E6| ) resulting in (^Sj), that is identical to the prefactor of the linear term in r 
in ( |D4]) (see also (|102|) ). The term resulting from \S\'^ is 
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{{\S{^.^)?)) 



'dip {\D{(p, (p) + D{tt - ip,n - ip) + D{-n + (p,-TT + (p) + D{-(p, -(p)\' 
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+ \D{<p, -<p) + D{n - ip, -TT + <p) + D{-<p, <p) + L)(-7r + ^, vr - (^)p 
+ \D{ip, 7i-ip) + D{-ip, -7r + ip) + D{-7r + ip, -ip) + D{-k - ip, ip)]"^ 
+ \D{ip, -'K + ip) + D{-ip, Tc -ip) + D{-7i + ip,ip) + D{tx - ip, -</?)P) 

and it is identical to the prefactor of in (|D4|) . The last term is given by 



(E14) 



--{{m\ip,ip))) 



1 



327r Jo 



dip {[D{ip, ip) + D{'n-ip,'n-ip) + D{-Ti + ip, --n + ip) + D{-ip, -</?)]' 
+ [D{ip, -ip) + D{n - ip, -n + ip) + D{-ip, ip) + D^-n + (/?, vr - ip)]^ 
+ [D{ip, TT - ip) + D{-ip, -7f + ip) + D{-7i + ip, -ip) + D^-K - ip, ip)f 



+ [D{ip, -Tc + ip) + D{-'K + ip,ip) + D{-ip, TT - ip) + D{tt - ip, -ip)y 



(E15) 
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It is straightforward to show that ([E15| ) is identical to the prefactor of in (pil 



We have demonstrated that the results of |^ are in agreement with our results for all the cases 
which we investigated. 
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